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SECTION  1 
BACKGROUND 

The  recognition  by  the  Air  Force  Avionics  Laboratory  (AFAL) 
of  the  need  for  improved  sensors  and  techniques  for  implementing 
fire  control  director  algorithms  has  resulted  largely  from  the 
conclusions  of  two  programs  - EXPO  V and  Fire/Fly.  EXPO  V 
is  the  latest  in  a series  of  man-in-the-loop  simulation  studies 
of  new  fire  control  and  gunsight  concepts.  The  Fire/Fly 
program  is  the  integrated  Fire  Control/Flight  Control  study  ^ 

which  is  investigating  methods  for  integrating  the  fire  control 
system  with  the  flight  control  system  to  improve  effectiveness 
while  increasing  survivability.  Both  programs  have  cited  the 
necessity  for  improved  sensors  as  well  as  improved  director 
algorithms . 

The  purpose  of  the  Fire  Control  Systems  analysis  effort, 
summarized  in  the  following  report,  was  to  examine  existing  fire 
control  systems  test  data  (Sight  Eval) , identify  error  sources, 
evaluate  and  modify  fire  control  algorithms,  and  perform 
programming  for  simulation,  weapon  system  investigation,  and 
validation. 

The  tasks  contained  in  this  volume  were  of  an  analytical 
nature . 


1 


SECTION  2 
TASK  DEFINITION 


The  broad  categories  of  effort  contained  in  the  Statement 
of  Work  were  refined  by  coordination  with  Captain  J.  Silverthorn, 
RWT-2,  and  resulted  in  the  following  series  of  tasks. 


2.1  IDENTIFY  RADAR  RANGE  ERRORS  AND  DEVELOP  FILTER 

An  important  input  to  the  director  gunsight  computations 
is  target  range.  As  a result,  it  is  desirable  to  have  as 
accurate  a range  measurt^ment  as  possible.  Some  type  of  filter 
(Kalman  or  other)  is  needed  to  produce  an  accurate  estimate  of 
range  (and  range  rate) . To  accomplish  this  task,  an  approximate 
noise  model  is  needed  to  aid  in  the  filter  design  process. 

This  noise  model  can  also  be  used  in  the  Air  Force  Flight 
Dynamics  Laboratory  (AFFDL)  LAMARS  simulation. 

2.2  TARGET  STATE  MEASUREMENT 

Two  important  subsystems  of  the  director  system  are  the 
target  position  sensor  and  Kalman  angle  filter.  An  appropriate 
performance  measure  for  these  subsytems  is  their  ability  to 
estimate  target  state  (position,  velocity,  acceleration). 

As  a result,  to  evaluate  their  performance  it  is  necessary 
to  independently  measure  these  target  parameters  more  accurately 
than  the  expected  estimation  errors.  A technique  must  be 
developed  and  analyzed  to  accomplish  this  task. 

2.3  IDENTIFY  AND  CORRECT  RADAR  LAG  1PR0BLEM 

The  MA-1  radar  was  seen  to  lag  the  target  by  as  much  as 
50  mrad  during  the  Sight  Eval  flight  test.  Since  this  same 
radar  is  to  be  used  during  the  director  flight  test  for  initial 
acquisition  and  handoff  to  the  ASCOT  sensor,  it  is  necessary 
to  be  able  to  predict  the  .]Lag  so  that  the  net  error  (difference 
between  target  true  position  and  radar  estimate)  is  less  than 
ASCOT  scan  size  (+  15  mrad) . 
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SECTION  3 
ANALYSIS  SUPPORT 


3.1  IDENTIFY  RADAR  RANGE  ERRORS  AND  DEVELOP  FILTER 

3.1.1  Introduction 

This  report  evaluates  the  performance  of  certain  optimal 
and  sub-optimal  Kalman  filters  in  yielding  estimates  of  target 
state  variables.  The  evaluation  uses  simulated  trajectories 
with  noise  added  to  the  true  range  as  input  to  the  filters,  and 
compares  filter  estimates  of  the  state  variables  to  their  true 
values . 

A Kalman  filter  is  often  designed  with  respect  to  the 
coordinate  system  in  which  the  target  position  is  measured;  nor- 
mally the  inertial  or  line-of-sight  coordinate  system.  The  filters 
considered  in  this  report  use  only  range  information.  Therefore, 
to  give  a complete  specification  of  the  target  state,  the 
relative  range,  relative  radial  velocity,  and  relative  radial 
acceleration  estimates  yielded  by  these  filters  must  be  converted 
to  a useful  frame  by  coordinate  transformations  based  on  angular 
position,  angular  rate,  and  angular  accelerations  as  measured 
by  suitable  devices  and  estimated  by  other  filters. 

Four  Kalman  range  filters  were  compared.  The  first  three 
were  based  on  a constant  relative  radial  acceleration  model 
and  the  fourth  on  a constant  relative  radial  velocity  model. 

To  evaluate  and  compare  these  filters,  trajectories  of 
various  types  were  generated,  noise  was  added  to  each  true  range, 
and  the  resulting  "measurement"  was  used  as  input  to  each  filter. 
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Results  are  given  graphically  and  as  root-mean-square  deviations 

of  the  filter  estimate  of  range,  (radial)  velocity,  and  (radial) 

acceleration  over  a chosen  interval. 

The  simulated  constant  acceleration  trajectories  assume 

_2 

(constant)  accelerations  of  0,  Ig  (9.8  m s ) , and  6g  (58.8 
-2 

ms  ) . Two  simulations  assume  a constant  negative  acceleration. 
The  initial  range  and  velocity  were  assumed  to  be  600  m and  200 
m s respectively. 

Two  similar  trajectories  resulting  from  variable  acceler- 
ation were  also  used  as  simulations.  These  resulted  from 
sinusoidal  accelerations  of  8 sec.  period  and  amplitudes  of 
either  2g  or  6g.  The  initial  range  was  assumed  to  be  600  m and 
the  initial  velocity  was  zero. 


3.1.2  Theory 

Filter  A;  Constant  Acceleration  Model  - Time  Varying  Gain 


Let  the  state  vector  representing  the  target  be  denoted  by 

(t) 

x(t)  = 

x^  (t) 

where  x^(t),  X2(t),  and  represent  the  range,  (radial) 

velocity,  and  (radial)  acceleration,  respectively,  of  the  target 
at  time  t.  The  constant  acceleration  model  assumes  that  the 
state  vector  evolves  in  time  according  to 

Xi=X2 

X2=X3 

x^=w(t) 

where  w(t)  is  assumed  to  be  zero-mean,  Gaussian  noise  uncorrelated 
with  itself  in  time.  This  model  can  be  written  in  matrix  form  as 


X(t)  = Fx(t)  + Gw(t) 


(1) 


where 


o 

O 

t 

’ 0* 

F = 

0 0 1 

. G = 

0 

0 0 0 

1 

^ * 

. 

Assuming  that  w(t)  is  constant  over  the  interval  (to/t), 
the  solution  of  Eq.  (1)  is  given  by 


x(t)  = <<>  (t,  t^)x(t^)  (t,tT)G(T)dT  jw{t) 


where 


<l>(t,T)  = 


1 

0 

0 


1 2 
t-T  i(t-T)^ 


(2) 


1 

0 


t-T 

1 


For  successive  times  and  t^.  Eg.  (2)  becomes 

.t. 


^^k+1'  ^k^-^^k^ 


J. 


'k+l 


w{t,^)  (3) 


Defining 


x(tk^l) = x(k+l) 
x(tj^)  = x(k) 
w(tj^)  E w(k) 


/ 


'k+l 


<|)  (tj^+l,T)G(T)dT  = r(k+l,k), 


Equation  (3)  can  be  written  as 

x(k+l)=(t)  (k+l,k)x(k)  + r(k+l,k)w(k) 
By  direct  integration. 
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Both  T and  <}>  are  constant  given  a fixed  interval  between 
successive  measurements. 


It  is  assumed  that  the  maneuver  noise  term,  w(k) , is 
zero-mean  Guassian  noise  with  the  autocorrelation  property 

E{w(j)w'  (k)}  = Q(k)6jj^  (4) 

where 

2 

Q(k)=a  (maneuver-noise  variance) 

w 

and  6 is  the  Kronecker  delta. 

A measurement  model 

z(k+l)  = Hx(k+1)  + v(k+l)  (5) 

is  assumed  where  v(k+l)  is  the  (scalar)  measurement  noise, 

H=[l  0 0] 

is  the  measurement  matrix,  and  z(k+l)  is  the  measurement  at 
time  Equation  (5)  may  be  written  in  the  scalar  form 

z(k+l)  = Xj (k+1)  + v(k+l)  (6) 

It  is  assumed  that  the  measurement  noise  has  the  properties 
E|v(j)v(k)}  = R(k)6j^  (7) 

and 

E|v(j)w(k)}  = 0 
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where 


R(k)  = aj 

is  the  measurement  noise  variance;  also  called  the  signal  noise 
variance . 

The  "best"  estimate  of  the  state  vector  at  time  tj^  is 
denoted  by  x(klk);  this  estimate  based  on  k measurements  z(k),  z(k-l), 
z(2),  z(l).  The  covariance  matrix  P(j|k)  is  defined  to  be 

P(j lk)=E  \ [x( j)-x( j I j) ] [x(k)-x(k|k) ] •[  , 

where  the  transpose  is  indicated  by  a prime. 

ViTith  these  definitions,  the  Kalman  equations  are,  in  the 

logical  order  of  their  execution 

P(k+l|  k)=(j)P(k|k)(t)'+rQ(k)r'  (i)  I 

K(k+l)=P(k+llk)H'  tHP(k+l|k)H'+o^]”^  (ii)  ( 

, ) (8) 

x(k+llk+l)=(l)X(k|k)  + K(k+1)  t z(k+l)-H(J)x(k|k)  ] (iii)y 
P(k+llk+l)=[I-K(k+l)H)p(k+lik)  (iv)  ) 

In  order  to  use  Eqs .( 8) , initial  estimates  for  S(k|k)  and 
P(k|k)  for  some  value  of  k are  required. 

If  the  position  of  the  target  is  known  exactly  at  three 
successive  times  tj^,  k = 1,  2,  3,  and  the  target  was  accelerating 
uniformly,  the  velocity  X2(3)  could  be  found  from  the  relation 

X2(k)  = ^ [3x^(3)  - 4x^(2)  + x^(l)]  (9) 

which  may  be  derived  from  elementary  kinematics.  If 
instead  the  measured  values  of  the  range  in  Eq.  (9)  are  used 
an  estimate  of  the  velocity  at  k = 3 is  given  by 
X2(3|3)  = ^ [3z(3)-4z(2)  + z (1)1 
Similarly,  the  acceleration  x^(3)  would  be  given  by  (again  from 
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elementary  kinematics) 


X3(3|3)  = [z(3)  - 2z(2)  + z(l)] 

The  first  value  of  k for  which  these  estimates  may  be 
formed  is  k = 3;  therefore,  the  earliest  estimate  of  the  state 
vector  which  can  be  made  using  the  above  procedure  is 


X3(3|3)  = 


[3z(3)-4z(2)  + z(l)  ] 


[z(3)-2z(2)  + z{1)  ] 


Because  of  the  extreme  sensitivity  of  X3(3|3)  to  noise, 
better  (and  simpler)  estimate  is  given  by 


x(3|3)  = 


[3z (3)-4z (2)+z (1) ] 


To  form  an  initial  estimate  of  P(3|3),  we  will  use  Eq. 
(10) . The  state  vector  at  k = 3 is  represented  by 


Xj^  (3) 


x(3)  = 


^ [3x^(3)-4x3(2)+x^(1)] 


[x^(3)-2xj^(2)+xj^(l)  ] 


Combining  Eqs.  (6),  (10),  and  (12), 


X3  (3) 


1 


x(3)-x(3|3)  = ^ [3x3(3)-4x^(2)+Xt  (1)  ] [ 3 z ( 3) -4  z ( 2 ) +z  ( 1)  ] 


[x,  (3)-2x,  (2)+x,  (1)  ] 


[z(3)-2z(2)  + z(l)  ] 


v(3) 


^ [3v(3)-4v(2)+v(l) ] 
1 


[v(3)-2v{2)+v(l) ] 


From  its  definition  and  from  Eqs.  (4)  and  (7),  P(3l3) 
is  given  by 


P(3| 3)=E 


9 

p- 
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1 
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12T 

.2T^ 

12T 

12  . 
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VJith  Eqs.  (11)  and  (13)  as  initial  estimates  of  the  state 
vector  aii-d  covariance  matrix,  the  state  vector  estimate  is 
updated  with  each  new  measurement  using  the  Kalman  Eqs.  (8).  [1]. 


Filter  B;  Constant  acceleration  model.  Synthetic  time-varying  gain. 
The  Kalman  gain,  K(k) , computed  from  the  second  equation 


of  the  set  (8)  of  the  previous  section,  is  a function  of  the 


measurement  interval  T,  the  measurement  noise  variance  a , the 

V 


maneuver  noise  variance  a^,  and  the  iteration  number  k.  To  apply 


the  Kalman  equations  to  a particular  filtering  problem,  these 
parameters  must  be  known,  and  the  Kalman  gains  may  be  calculated 
for  each  k independent  of  any  measurement.  The  components  of 
K(k)  are  therefore  determined  a priori. 

A measurement  interval  of  T = 0.04  seconds  has  been  used 
throughout  all  simulations.  For  Filter  A values  of  K(k)  from 


11 


r 

i 

j 

! 


k = 3 to  k = 200  were  computed.  It  was  supposed  that  each  com- 
ponent of  the  Kalman  gain  could  be  approximated  by  a function  of 
the  form 

' ^i 

K (k)  — + K.  (oo)  , i=l,2,3  (14) 

k -b.  ^ 

1 


where  the  parameters  a^  and  b^  are  found  from  fitting  Eg.  (14) 
to  the  values  K^(5)  and  K^(30)  which  were  computed  for  Filter 
A.  K^(oo)  is  the  value  to  which  K^(k)is  an  asymptote.  Although 
Eqs.  (14)  do  not  purport  to  be  the  best  possible  fit  to  the 
optimum  gain  components,  the  advantage  is  that  the  (k)  are 
easy  to  calculate. 

If  Eg.  (14)  is  used  to  calculate  the  Kalman  gain  components, 
the  state  estimate  is  updated  according  to  the  single  equation 
x(k  |k)=<|)x(k-l  |k-l)+K'  (k)  [ z (k) -H())X (k-1  |k-l)  ] 
with  a considerable  saving  in  computation  time. 


Filter  C;  Constant  Acceleration  Model.  Constant  gain.  1 

Figures  (1),  (2),  and  (3)  show  that  each  component  of  the 

Kalman  gain  assumes  its  largest  value  at  k = 4,  and  decreases 
monotonically , asymptotically  approaching  some  minimum  value. 

The  motivation  of  this  filter  is  to  find  some  particular  value 
of  k,  say  k = i,  such  that  the  constant  gain  K(i)  may  be  used 
to  update  the  state  estimate  for  each  iteration;  i.e., 

X (k+1 1 k+1)  =())x  (k  I k)  +K  (i)  [ z(k+l)  -H(()x  (k  | k)  ] 

If  i is  chosen  to  be  small,  the  Kalman  gain  will  be 
relatively  large  and  will  remain  so,  since  K(i)  is  constant, 
by  hypothesis.  The  resulting  estimates  will  be  very  noisy  and 


show  little  tendency  to  converge.  If  the  value  of  i chosen  is 
large,  the  corresponding  Kalman  gain  will  be  small  and  errors 
in  the  initial  estmiates  will  be  slowly  corrected.  It  is  hoped 
that  an  intermediate  value  of  i can  be  found  such  that  after  a 
reasonable  number  of  iterations,  convergence  can  be  achieved. 

Filter  D:  Constant  Velocity  Model.  Time  varying  gain. 

As  its  name  suggests,  this  filter  assumes  that  the  target 
is  moving  at  constant  velocity.  The  development  of  the  filter 
equations  for  this  model  is  parallel  to  that  of  Filter  A.  It 
is  simpler  in  the  sense  that  the  state  vector  has  two  components 
rather  than  three.  The  associated  matrices  are  hence  2x2, 
rather  than  3x3. 

Let  the  state  vector  of  the  target  be  denoted  by 


x(t)  = 


(t) 

X2(t) 


where  Xj^(t)  represents  the  range  of  the  target  and  X2(t)  is  its 
(radial)  velocity.  This  model  assumes  that  the  state  vector 
evolves  according  to 


Xi  = X2 

X2  = w(t) 


where  w(t)  is  assumed  to  be  zero-mean,  Gaussian  noise  uncorrelated 
with  itself  in  time.  In  matrix  notation, 

x(t)  = Fx(t)  + Gw(t)  (15) 


W 


where 


o 

1 

0 

F = 

, G = 

I 

O 

o 

1 

1 

The  solution  of  Eq.  (15)  is  given  by 
x(t)  = (J)(t,tQ)x(t^)  + (})(t,  t)G(  T)dT]w(t) 


where 


<{)(t,T)  = 


o 

1 t-  t) 


becomes 


0 1 

For  successive  measurement  times  and  tj^,  Eq.  (16) 

^k+1 


/K-l-J. 

())(tK+i,T  )G(T)dT]  W(tj^) 


Defining,  as  before 
x(tk^l) = x(k+l) 


x'tj^)  = x(k) 

/^k+1 

-k 


^H+1'  ")G('r)dTEr(k+i;  k) 

VJith  the  above  definitions,  Eq.  (16)  may  be  written  as 


X (k+1)  =4)  (k+1 , k)x(k)  + r(k+l,  k)w(k) 
and  by  direct  integration, 

r(k+l,  k)= 

where 


^ktl  - t 


r^k+1 

1 m2' 

j 't’^^k+l'  T)G(T)dT  = 

— T 

2 

ft 

T 

k 

(16) 


(17) 
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and 


(j)  (k+1,  k)  “ (t)= 


which  is  constant,  given  a fixed  interval  between  successive 
measurements . 


Again,  assume  that 


E I w(  j)  w(k)  } = Q(k)  6 . 


where 


Q(k)  =a; 


(maneuver  noise  variance) 


The  measurement  model  is,  as  before 


z(k+l)  = Hx (k+1) +v (k+1) 

where  v(k+l)  is  the  (scalar)  measurement  noise  and 
H = [ 1 0 ] 

so  that  eq.  (18)  may  be  written 

z(k+l)  = Xj^  (k+1) +v  (k+1) 


Again,  assume  that  the  measurement  noise  has  the  properties 


E I v(  j)  v(k)  } =R(k) 


R(k)=a; 


E I v(^j  ) w(k)}  =0 

I 

The  Kalman  equations  for  this  model  have  the  same  +orm  as 
Eq.  (8) , where  all  quantities  in  those  equations  assume  the 
definitions  given  in  this  section.  The  initial  state  estimate 
and  covariance  matrix  are  initialized  for  k = 2 as  follows: 


r 


and 


x(2  2)  = 


P(2|2)  = ^2 


z(2) 

i [z(2)-z(l)] 


If  :1 


(19) 


(20) 


P(2|2)  is  obtained  in  a manner  analogous  to  the  procedure 
used  for  Filter  A. 

VJith  these  initial  estimates  of  the  state  vector  and 
covariance  matrix,  the  state  vector  estimate  is  updated  with 
each  new  measurement  by  means  of  the  Kalman  Eqs.  (8).  [1]. 


3.1.3  Performance 

Two  types  of  trajectories  were  used  in  the  simulations: 

those  resulting  from  a constant  acceleration  (called  constant 

acceleration  trajectories)  and  those  resulting  from  a sinusoidal 

acceleration  (called  sinusoidal  acceleration  trajectories) . 

A constant  acceleration  trajectory  conforms  to  the  model  of 

Filters  A,  B,  and  C,  but  is  a departure  from  the  model  of  Filter 

D,  except  for  x^d^)  =0.  A sinusoidal  acceleration  trajectory 

is  a departure  from  the  models  of  all  filters.  This  type  of 

2 

trajectory  serves  to  determine  the  value  of  (maneuver  noise 
variance)  which  will  optimize  the  response  of  the  filter  in  the 
sense  that  it  will  track  departures  from  the  model,  yet  not  be 
too  sensitive  to  measurement  noise. 

A trajectory  is  completely  determined  by  the  form  of  the 
acceleration,  X2(k),  the  initial  range  x^(0),  and  the  initial 
velocity  X2(0).  Range,  velocity,  and  acceleration  are  then 


Ik 
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r 


computed  for  each  k.  After  generating  the  true  values  Xj^(k), 

X2  (k) , and  X2(k),  a noise  term  is  computed  as  follows. 

Twelve  random  numbers,  each  between  0 and  1,  are  generated 
and  added.  The  number  6 is  then  subtracted  from  the  sum.  The 
result,  when  multiplied  by  is  an  approximately  zero-mean, 

Guassian  random  variable  with  standard  deviation  a . In  symbols, 

( 12  \ 

v(k)=a^|  %"k 


where  N.  is  a random  number  between  0 and  1. 

1 

Starting  with  k = 1,  the  same  sequence  of  random  variables 
nj^  was  generated  for  all  runs.  All  variable  acceleration 
simulations  used  the  same  sequence  of  random  variables;  i.e., 
for  a given  k,  the  value  of  was  always  the  same.  The 
constant  acceleration  trajectories  use  the  first  50  values  of 
of  this  set  to  compute  the  noise  term,  v(k). 

Figure  1 is  a plot  of  v(k)  for  k = 1 to  50  for  a standard 
deviation  ^^^=6  m.  Table  I lists  the  values  of  v(k)  from  k = 1 
to  50  for  0^=6  m.  Note  the  preponderance  of  negative  noise  terms 
in  the  last  few  entries.  The  average  of  the  last  eight  noise 
terms  is  -2.59  m and  the  average  of  the  last  six  is  -2.99  This 

results  in  a negative  bias  in  the  estimates  for  the  last 
several  terms  of  the  constant  acceleration  simulations. 

Measurement  noise  characterized  by  cr^=ij  2,  4,  6,  8,  10  m[ 
was  used  in  these  simulations.  A study  by  Oricon  [2]  suggests 
a value  of  cr^=3  m for  radar  of  "medium  accuracy",  but  0^=6  m was 
chosen  as  a representative  of  the  measurement  noise  encountered 
at  range  of  about  600  m.  Section  3.1.3  of  this  report  shows 
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TABLE  1 


NOISE  TERMS  v(k) 
o =6m 

V 


k 


v(k) 


1 

-14.75 

2 

- 6.92 

3 

3.74 

4 

5.56 

5 

11.82 

6 

- 6.24 

7 

4.24 

8 

4.65 

9 

.40 

10 

15.38 

11 

- 6.57 

12 

- . 99 

13 

2.20 

14 

8.97 

15 

- 5.79 

16 

- 1.33 

17 

5.55 

18 

- 5.28 

19 

-10.98 

20 

- 1.42 

21 

1.87 

22 

1.12 

23 

3.37 

24 

1.08 

25 

2.88 

26 

- 7.73 

27 

- .15 

28 

8.57 

29 

9.49 

30 

5.68 

31 

- 3.06 

32 

5.00 

33 

- 1.37 

34 

2.81 

35 

11.02 

36 

4.19 

37 

.53 

38 

3.00 

39 

- 4.03 

40 

15.76 

41 

8.69 

42 

2.81 

43 

- 3.33 

44 

.70 

45 

- 6.93 

46 

- 6.49 

47 

- 5.18 

48 

3.20 

49 

4.66 

50 

19  - 7.19 

region  of 
negative 


bias 


to  be  about  6 m for  the  MA-1  radar  and  to  be  only  weakly 
dependent , 

Table  2 is  a summarv  of  simulated  traiectories  and  noise 
variances  used  to  evaluate  the  filters.  The  line  numbers  at 
the  left  identifv  the  traiectory  and  the  filter  parameters  used. 

At  the  right  of  the  table,  the  letter  A,  B,  C,  or  D indicates 
the  filters  used  to  estimate  state  variables  from  the  specified 
trajectory.  A hyphen  (-)  indicates  the  filter  was  not  used, 
for  reasons  which  will  be  presented  later. 

Plots  of  filter  performance  are  to  be  found  in  Appendix 
A.  Each  figure  in  this  appendix  is  a plot  of  the  true  velocity 
and  acceleration  and  their  estimates  by  each  of  the  four  filters. 
The  true  velocity  and  acceleration  are  plotted  with  solid  lines, 
and  their  corresponding  estimates  are  plotted  with  dashed  lines. 
Since  Filter  D does  not  yield  an  acceleration  estimate,  no  such 
plot  will,  of  course,  be  present.  Since  Filter  C was  used  only 
in  cases  16,  17,  18,  24,  and  30,  estimates  from  this  filter  will 
be  present  in  these  cases  and  absent  in  all  others. 

The  format  of  each  figure  is  as  follows.  The  numerical 
part  of  the  figure  number  is  the  same  as  the  case  number  in  Table 
2.  The  plots  in  the  left  column  of  each  figure  show  X2(k)  and 
X2(k|k)  vs  k for  each  filter  used.  There  will  be  either  three 
or  four  such  plots  in  this  column,  depending  on  whether  or  not 
Filter  C was  used  for  that  case.  The  right  column  of  each  figure 
shows  plots  of  x^Ck)  and  x^Ckjk)  vs  k for  Filters  A and  B and  for 
Filter  C as  well  in  those  cases  where  this  filter  was  used.  There 
will  be  either  two  or  three  such  plots  in  the  right  column. 
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Case 


TABLE  2 


Range 


Simulated  Trajectories  and  Noise  Vaxiaiicea 


Velocity 


Acceleration 


Measurement 
Noise  Variance 


Maneuver 

Noise 

Variance 


Filter 


X,  (0) 

^ (m) 

X 

2 -1 
m s 

(k) 

-2 

m s 

0 2 

V 

0 2 
w 

1 

600 

200 

0 

4 

0 

A 

B 

- 

D 

2 

II 

It 

9.8 

It 

It 

A 

B 

- 

D 

3 

II 

II 

58.8 

It 

It 

A 

B 

- 

D 

4 

tl 

tl 

0 

16 

tl 

A 

B 

- 

D 

5 

II 

II 

9.8 

II 

tt 

A 

B 

- 

D 

6 

II 

It 

58.8 

tl 

tt 

A 

B 

- 

D 

7 

tl 

It 

0 

36 

tl 

A 

B 

- 

D 

8 

tl 

It 

9.8 

tl 

tl 

A 

B 

- 

D 

9 

tl 

It 

58.8 

It 

tt 

A 

B 

- 

D 

10 

II 

It 

0 

64 

tt 

A 

B 

- 

D 

11 

It 

It 

9.8 

tl 

ft 

A 

B 

- 

D 

12 

tl 

It 

58.8 

II 

It 

A 

B 

- 

D 

13 

It 

It 

0 

100 

tl 

A 

B 

- 

D 

14 

tl 

It 

9.8 

tl 

It 

A 

B 

- 

D 

15 

It 

tt 

58.8 

tl 

It 

A 

B 

- 

D 

16 

It 

II 

0 

36 

6000 

A 

B 

c 

D 

17 

rr 

tl 

9.8 

tt 

It 

A 

B 

C 

D 

18 

tl 

It 

58,0 

It 

tt 

A 

B 

c 

D 

19 

It 

II 

-9.8 

1! 

0 

A 

B 

- 

D 

20 

tl 

-58.8 

tl 

It 

A 

B 

- 

D 

21 

If 

II 

36 

0 

A 

B 

- 

D 

22 

It 

tl 

tl 

II 

2000 

A 

B 

- 

D 

23 

tl 

tl 

tl 

tt 

4000 

A 

B 

- 

D 

24 

It 

It 

It 

tl 

6000 

A 

B 

c 

D 

25 

It 

tl 

II 

tl 

8000 

A 

B 

- 

D 

26 

tl 

It 

II 

. TT  Ic 

It 

10000 

A 

B 

- 

D 

27 

II 

It 

^""-100 

tl 

0 

A 

B 

- 

D 

28 

II 

It 

tl 

It 

2000 

A 

B 

- 

D 

29 

tl 

tt 

II 

It 

4000 

A 

B 

- 

D 

30 

It 

II 

It 

It 

6000 

A 

B 

c 

D 

31 

It 

tt 

ft 

tt 

8000 

A 

B 

- 

D 

32 

It 

It 

tt 

tt 

10000 

A 

B 

— 

D 

21 


1 

1 


1 

I 
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Thus,  Figure  A-15  contains  plots  of  the  performance  of  Filters 
A,  B,  and  C for  the  trajectory  (and  noise  variances)  listed  under 

Case  15;  i.e.,  Xj^(O)  = 600  m,  X2(0)  = 200  m,  x^Ck)  = 58.8  m s 

2 2 2 
a = 100  m , a = 0. 

V w 

Cases  1 through  20  of  Table  2 list  constant  acceleration 

simulations.  These  conform  to  the  model  of  Filters  A,  B,  and  C. 

Cases  1,  4,  7,  10,  13,  and  16  also  conform  to  the  model  of  Filter  D. 

2 

In  Cases  1 through  15,  a =0  and  measurement  noise 

w 

variances  are  given  the  values  = {4,  16,  36,  64,  100  m^  } . 

2 2 2 2 “ 6 

In  Cases  16  through  18,  a = 36  m and  cr  = 6000  m s ; the 

^ V w 

values  judged  "best"  to  handle  deviations  from  the  model.  Data 
supporting  this  conclusion  are  listed  in  Section  3.1.4.  Cases 
19  through  20  specify  negative  accelerations.  Cases  21  through  26 
specify  sinusoidal  accelerations  of  2g  amplitude  and  8 s period. 

Cases  27  through  32  specify  sinusoidal  accelerations  of  the  same  ' 

period,  but  with  a 6g  amplitude. 

Filter  A;  Constant  Acceleration  Model.  Time  varying  gain. 

Figure  A-1  through  A-20  contain  plots  of  the  velocity  and 

acceleration  estimates  of  Filter  A for  Cases  1 through  20  listed 

2 

in  Table  2.  The  value  of  was  used  in  both  the  noise  aenerator  ; 

and  in  the  filter.  These  plots  show  both  velocity  and  acceleration 
estimates  to  be  extremely  noisy  for  the  first  half-second,  but 

! 

they  smooth  out  considerably  at  the  end  of  one  second,  and  j 

approach  the  true  values  quite  closely  after  about  30  iterations  ! 

(1.2  s) . As  expected,  the  lower  the  value  of  a^,  the  faster  and  j 

closer  the  convergence.  Note  the  characteristic  "tail",  or 


22 


negative  bias,  in  the  estimates,  particularly  those  of  higher 
signal  noise  variance.  This  is  caused  by  the  anomoly  in  the 
random  numbers  referred  to  earlier. 

Figure  A-21  through  A-32  contain  plots  of  velocity  and 
acceleration  estimates  of  sinusoidal  acceleration  simulations. 

Several  values  of  the  maneuver  noise  variance  are  used  to  find 

2 

the  "best"  value  of  a . (See  Section  3.1.4.)  It  is  clear  from 

w 

. 2 

Figure  A-21  and  A-27  that  filter  estimates  using  = 0 are 

extremely  poor,  as  one  would  expect. 

Filter  B;  Constant  Acceleration  Model.  Synthetic  gain. 

The  motivation  of  this  filter  and  the  method  of  calculation 

of  the  parameters  used  to  generate  K' (k)  are  discussed  in  Section 

3.1.2.  Since  the  values  of  K(k),  the  optimum  time-varying  gain, 

2 2 

depend  upon  the  noise  parameters  and  a^,  the  parameters  which 
generate  K'  (k)  were  found  by  fitting  Eq.  (14)  to  K^(5)  and  Kj^(30) 
for  each  pair  of  noise  parameters. 

I 

Figures  2,  3,  and  4 are  plots  of  K^(k)  and  K^(k)  vs  k for 

the  three  components  of  the  Kalman  gain.  The  noise  parameters 

used  to  generate  K(k)  were  those  judged  to  be  optimum  in  a sense 

2 2 2 2 “ 6 

discussed  in  the  next  section;  a = 36  m and  a = 6000  m s 

V w 

Although  the  fit  is  not  a close  one,  this  procedure  yields  synthetic 
Kalman  gains  such  that  the  estimates  of  Filter  B in  Figures  A-1 
through  A-32  are  difficult  to  distinguish  from  corresponding  plots 
of  Filter  A. 

Filter  C;  Constant  Acceleration  Model.  Constant  gain. 

Figures  2,  3,  and  4 are  plots  of  the  components  of  the 
optimum  Kalman  gain  vs  the  iteration  number  k.  The  optimum  constant 
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Figure  2.  Position  Kalman  Gain  Component 


1 

j 

value  of  the  Kalman  gain  K(i)  was  determined  by  finding  the  gain  | 

which  minimizes  the  rms  deviation  of  the  velocity.  (See  Section  j 

1 

3.1.4.)  This  gain  was  found  to  be  K(24),  and  is  shown  in  Figures 
2,  3,  and  4 as  a dotted  horizontal  line.  Since  these  values  were 
based  on  optimum  gains  for  - 36  m and  = 6000  m'^s’*  , we 
used  Filter  C only  where  appropriate;  i.e.,  for  parameters  listed 
on  lines  16,  17,  18,  24,  and  30.  Figure  A-16,  A-17,  and  A-18 
show  that  velocity  estimates  of  Filter  C tend  to  converge  to  the 
true  values,  though  not  as  quickly  nor  as  closely  as  the  corresponding 
velocity  estimates  of  Filters  A and  B.  Figure  A-16,  A-17,  and 
A-18  show  that  the  acceleration  estimates  of  Filter  C also  show 
a tendency  to  converge,  but  more  weakly  than  those  of  Filters  A 
and  B.  Figure  A-24  and  A-30  show  noisier  velocity  estimates  than 
those  shown  in  Figures  A and  B,  and  the  acceleration  estimates 
shown  in  these  figures  are  much  noisier. 

Filter  D;  Constant  Velocity  Model. 

This  filter  assumes  the  target  to  be  moving  with  constant 

relative  velocity.  Departures  from  constant  velocity  can  be 

2 

treated  by  choosing  a suitable  value  of  a^.  This  filter  yields 
position  estimates  and  velocity  estimates  only.  Attempts  to 
estimate  acceleration  based  on  some  algorithm  using  past  velocity 
estimates  proved  fruitless. 

It  can  be  seen  from  the  figures  that  even  for  very  large 
2 

values  of  a^,  this  filter  yields  a biased  velocity  estimate  even 
for  accelerations  of  only  Ig,  which  becomes  considerably  larger 
for  6g  accelerations.  The  velocity  estimates  are  quite  poor  for 
sinusoidal  acceleration  simulations,  though  better  than  Filter  C. 
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3.1.4 


Evaluation 


2 

Given  a trajectory  and  a pair  of  noise  variances  (a^  and 

2 

) we  wish  to  characterize  the  performance  of  each  filter  by  a 
single  number.  A root-mean-square  deviation  is  an  attractive 
candidate  for  such  a characteristic  number.  The  first  difficulty 
with  this  idea,  however,  is  that  the  filter  output  consists  of 
three  time-varying  quantities  x^(k.|k),  X2(k|k),  and  X2(k|k);  i.e., 
a range  estimate,  a velocity  estimate,  and  an  acceleration  estimate. 
The  smallest  rms  deviation  for  one  component  of  the  state  estimate 
may  not  be  the  smallest  for  the  others.  Another  difficulty  can  be 
seen  by  comparing  plots  of  the  estimates  of  Filter  A with  those 
of  Filter  D.  Although  Filter  A produces  noisy  velocity  and 
acceleration  estimates  for  low  k,  the  estimates  are  quite  good 
for  large  k (constant  acceleration  trajectories) , with  the 
exception  of  the  "tail"  referred  to  earlier.  Although  Filter  D 
yields  a less  noisy  velocity  estimate  for  low  k,  it  is  biased 
and  slowly  diverges  from  the  true  velocity  for  large  k.  Thus 
the  choice  of  the  better  filter  depends  upon  the  allowable  time 
between  acquisition  and  when  estimates  are  needed.  For  the  constant 
acceleration  trajectories,  we  have  assumed  this  interval  to  be 
1.5  sec.  Thus,  if  the  velocity  rms  deviation  were  computed  over 
the  2 sec.  interval,  one  would  find  the  deviation  lower  for 
Filter  D than  for  Filter  A,  due  to  the  high  noise  in  Filter  A for 
small  k.  Such  a procedure  would  lead  to  Filter  D being  judged 
"best".  But  if  one  has  at  least  1.5  sec.  between  acquisition 
and  the  time  an  estimate  is  needed,  a comparison  of  the  estimates 
of  Filter  A with  those  of  Filter  D shows  that  Filter  A is  best. 

A A 

since  X2(k|k)  and  x^ikjk)  of  Filter  A smooth  out  in  about  1.5  s. 
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it  v/as  decided  to  compute  the  rms  estimates  from  k = 39  (1.56  s) 
to  k = 50  (2,00  s) . In  symbols. 


(Ax.  ) 

1 rms 


[x^(k)-x^(klk) ] 


i=l,2,3 


where  the  subscript  i refers  to  one  of  the  three  components  of 
the  state  vector  and  the  state  estimate.  Results  are  listed  in 
Tables  3,  4,  and  5. 

For  sinusoidal  acceleration  simulations,  which  were  8 sec. 
in  duration  or  from  k = 1 to  k = 200,  the  figures  show  that 
Filters  A,  B,  and  D give  reasonable  estimates  after  about  3 sec. , 
so  the  rms  deviations  were  computed  from  k = 75  (3  sec.)  to 
k = 200  (8s);  i.e. , 


( Ax^ 


[x^(k)-x^(k|k)  ] 


i=l,2,3. 


Range  Estimates.  (Ax,) 

This,  of  course,  is  the  most  important  of  the  three 
estimates,  but  a reasonable  estimate  is  not  difficult  to  obtain. 

Cases  1 through  20  of  Table  3 list  the  rms  range 
deviations  for  constant  acceleration  trajectories.  It  can  be 
seen  that  Filter  A characteristically  has  the  smallest 
deviations,  with  Filter  B a close  second.  The  deviations  of 
Filter  C are  about  25%  higher  than  Filter  A.  The  deviations  of 
Filter  D are  quite  erratic;  small  for  small  values  of  the 
acceleration  and  much  larger  for  large  acceleration.  This  is 
to  be  expected,  for  Filter  D is  based  on  a constant  velocity 
mode 1 . 
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TABLE  3 RMS  Range  Deviations 
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TABLE  4 RMS  Velocity  Deviations 


TABLE  5 RMS  Acceleration  Deviations 
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Lines  21  throuqh  32  of  Table  3 list  (Ax,)  for  the 

1 rms 

sinusoidal  acceleration  simulations.  Trajectories  of  6g 
amplitude  represent  a more  extreme  departure  from  the  constant 
acceleration  model  than  are  likely  to  be  encountered  in  practice, 
and  are  included  to  test  filter  response.  Trajectories  of  2g 
amplitude  are  more  realistic,  and  therefore  greater  weight  is 
given  to  these  cases. 

Filter  B has  values  of  (Ax.)  only  a few  per  cent  greater 

1 rms  j ^ 

than  Filter  A,  and  the  deviations  of  Filter  D are  even  less 

than  those  of  Filter  B.  Filter  C contains  the  implicit  parameters 
2 2 2 2 ” 6 

= 36  m and  = 6000  m s , and  was  used  only  for  trajectories 
listed  on  lines  24  and  30.  The  deviation  is  notably  poorer 
than  that  of  the  other  filters. 

Velocity  Estimates.  (Ax„) 

2 rms 

Velocity  estimates  are  more  sensitive  to  noise  than  range 

estimates,  and  (Ax„)  will  vary  more  widely  for  different 

2 rms 

filters  than  (Ax.)  . Lines  1 through  20  of  Table  4 list 

1 rms  ^ 

(Ax-)  of  the  velocity  estimates  x-(k  k).  For  nearly  all 
2 rms  2 

I cases,  deviations  of  Filters  A and  B are  more  nearly  the  same. 

In  fact,  for  most  cases,  the  deviations  of  Filter  B are  slightly 
‘ less  than  those  of  Filter  A,  but  this  result  is  most  likely 

I 

anomolous.  Filter  D yields  erratic  deviations;  small  for 

trajectories  which  match  its  model  and  much  larger  for  those 

which  do  not.  Filter  C yields  consistently  large  deviations. 

Lines  21  through  32  list  (Ax„)  for  the  sinusoidal 
^ 2 rms 

acceleration  simulations.  As  above.  Filters  A and  B differ 
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but  slightly  and  Filter  C has  much  larger  deviations,  but  the 
deviations  of  Filter  D are  surprisingly  small,  though  not  so 
small  as  those  of  Filter  B. 

Acceleration  Estimates . ^ 3 ^ rms 

Lines  1 through  20  of  Table  5 list  Filter  A, 

B,  and  C.  There  are  no  deviations  computed  for  Filter  D,  for 

this  filter  does  not  provide  an  acceleration  estimate.  The 

deviations  of  Filters  A and  B are  nearly  the  same,  while 

those  of  Filter  C are  about  three  times  greater. 

Lines  21  through  32  list  for  sinusoidal  acceleration 

^ 3 rms 

simulations.  For  the  2g  amplitude  trajectories.  Filters  A and 
B have  deviations  of  about  Ig  and  Filter  C has  deviations 
about  twice  as  great. 

Selection  of  Maneuver  Noise  Variance 

2 

As  noted  earlier,  the  maneuver  noise  variance  (a  ) is 

w 

needed  to  allow  the  filter  to  respond  to  deviations  from  the 

model  on  which  it  is  based.  Of  all  simulations,  only  the 

sinusoidal  trajectories  (Cases  21  through  32)  are  deviations 

2 

from  the  models  of  all  filters,  and  so  the  best  value  of  is 

that  which  minimizes  (Ax„)  and  (Ax-)  for  these  cases. 

2 rms  3 rms 

Since  the  2g  simulations  are  more  typical  of  an  encounter  than 

2 

the  6g  simulations,  the  choice  of  a value  of  will  be  based 

upon  the  rms  deviations  resulting  from  Cases  21  through  26  of 

Tables  4 and  5.  For  Filter  A,  the  minimum  rms  deviations  for 

2 2 "*6 

both  (Ax_)  and  (Ax.,)  occurs  for  a,,  = 8000  m s 
2 rms  3 rms  w 
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2 2 —6 

For  Filter  B,  the  minimum  (Ax_)  occurs  for  a = 6000  m s 

2 rms  w 

2 2 “ 6 

and  the  minimum  {Ax_)  occurs  for  a = 8000  m s . Since  the 

3 rms  w 

differences  in  the  rms  deviations  resulting  from  these  two 

2 2 “6 

choices  are  quite  small,  = 6000  m s was  chosen,  for  the 
use  of  the  smaller  value  of  the  maneuver  noise  variance  will 
result  in  a slightly  less  noisy  filter. 

3.1.5  Conclusion 

The  most  striking  feature  of  the  tables  of  rms  deviations 

is  that  in  nearly  all  cases  the  deviations  of  Filters  A and  B 

are  within  a few  percen'-.  of  one  another.  Filter  A computes  the 

Kalman  gain  components  from  Eqs.  (i)  and  (iv)  of  equation  set 

(8)  by  matrix  multiplications,  whereas  Filter  B computes  the 

gain  components  from  three  simple  algebraic  equations. 

2 2 

Filter  A requires  the  parameter  o^,  a^,  and  T in  order 
to  give  numerical  results.  If  these  parameters  are  given,  Eqs. 
(8i) , (8ii)  , and  (8iv) , together  with  Eq.  (13)  - initialization 

of  the  covariance  matrix  - can  be  used  to  calculate  K(k) , 
independently  of  any  measurements.  By  following  the  procedure 
described  in  Section  II,  the  parameters  a^ , b , and  K^(oo)  can 
be  found  for  i = 1,  2,  and  3.  The  Kalman  gains  can  be  quickly 
computed  from  Eqs.  (21)  and  these  nine  parameters  at  a fraction 
of  the  time  needed  by  the  iterative  procedure  of  Eqs.  (8)  to 
yield  estimates  of  only  slightly  lower  quality  than  those  of 
the  optimum  filter.  Filter  A. 

Although  the  calculations  employed  by  Filter  C are  even 
simpler  than  those  needed  for  Filter  B,  the  results  are  much 
poorer  and  it  is  doubtful  that  the  estimates  would  be  useful. 
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Filter  D cannot  produce  a useful  acceleration  estimate, 
and  its  performance  in  estimating  velocities  is  very  poor  at 
accelerations  higher  than  Ig. 

[1]  Stochastic  Optimal  Linear  Estimation  and  Control, 

J.S.  Meditch,  McGraw-Hill ( 1969)  . 

[2]  Advanced  Defensive  Fire  Control  Processing  Methods 
Study.  ORINCON  F33615-76-C-1303  Contract  Status 
Briefing  31  Jan  77. 
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3.2  TARGET  STATE  MEASUREMENT 

3.2.1  Introduction 

The  general  objective  of  the  target  state  measurement  task 
is  to  devise  independent  methods  for  evaluating  the  accuracy  of 
target  state  vectors  (position,  velocity,  and  acceleration) 
obtained  from  a director  fire  control  system.  The  particular 
director  system  to  be  considered  is  being  developed  by  Honeywell 
(the  integrating  contractor)  for  AFAL/RWT.  A simplified  block 
diagram  of  some  essential  features  of  this  system  is  shown  in 
Figure  5 ; components  and  subsystems  which  would  not  affect  an 
actual  combat  encounter  (e.g.,  the  Mini-Heads-Up  Display,  used  to 
obtain  film  records)  are  not  included  in  this  diagram.  Director 
fire  control  systems  are  characterized  by  an  extensive  use  of 
data  in  a line-of-sight  coordinate  frame  rather  than  an  aircraft 
body  frame,  and  this  characteristic  may  lead,  in  principle,  to 
a more  stable  gunsight  system.  Figures  shows  that  the  attacker's 
Inertial  Sensor  Assembly  and  Air  Data  Computer  (via  the  MA-1 
Fire  Control  System)  are  inputs  to  the  director  system  in  addition 
to  the  target  sensors  and  the  Time  Code  Generator  and  Mode 
Control  Unit.  The  final  output  of  the  director  system  is  presented 
on  the  Heads-Up  Display. 

The  task  objective  stated  above  has  several  possible 
interpretation  features.  First,  the  stage  in  the  director 
system  at  which  target  position,  velocity,  and  acceleration  are 
to  be  measured  is  not  specified.  Presumably,  in  order  to  evaluate 
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' the  performance  of  the  entire  system,  the  most  appropriate  stage 

would  be  as  close  as  possible  to  the  final  output  of  the  system, 

, i.e.,  the  Head-Up  Display.  On  the  other  hand,  it  may  be  possible 

to  identify  inadequacies  in  various  subsystems  if  target  state 
vectors  are  evaluated  at  appropriate  states  closer  to  the 
sensor  inputs. 

A second  interpretation  feature  concerns  the  coordinate 
frame  in  which  the  "true"  target  state  vectors  are  to  be  com- 
pared with  target  state  vectors  from  the  director  system.  There 
are  at  least  three  commonly  used  reference  frames:  (1)  an 
inertial  frame,  fixed  in  the  earth  or  in  a local  air  mass,  in 
which  the  coordinate  directions  are  north(x),  east(y),  and  down(2); 

(2)  an  aircraft  body  frame,  in  which  the  coordinate  directions 
are  forward  along  the  aircraft  axis(x)*,  out  the  right  wing(y), 
and  the  right-hand  cross  product  of  these  directions ( z) ; and  (3) 
a line-of-sight  frame,  in  which  the  coordinate  directions  are  along 
the  line-of-sight  to  the  target (x),  toward  the  front  of  the  air- 
craft (y),  and  the  right-hand  cross  product  of  these  directions (z) . 

The  director  system,  of  course,  always  considers  the  target  state 
vectors  relative  to  the  attacker,  and  thus  the  evaluation  methods 
to  be  devised  should  also  refer  to  the  attacker  as  the  origin.  In 
general,  if  the  comparison  is  to  be  made  at  a stage  close  to  the 
target  sensor  inputs  of  the  director  system,  then  a line-of-sight 
frame  is  probably  indicated,  and  if  this  comparison  is  to  be  made 
at  a stage  close  to  the  final  output,  then  the  attacker's  body 
^ _ frame  is  probably  appropriate.  If  only  the  magnitudes  of  the 

j * A "velocity"  frame,  in  which  the  x direction  is  along  the  total 

) velocity,  is  also  used. 
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target  state  vectors  are  to  be  compared,  then  the  orientation 
of  the  coordinate  frame  need  not  be  specified. 

A third  interpretation  feature  concerns  the  independence 
of  the  evaluation  methods.  Completely  independent  methods 
would  not  make  use  of  any  of  the  sensor  inputs  or  computational 
algorithms  used  by  the  director  system.  However,  such  methods, 
although  possible  in  principle,  would  generally  be  difficult  or 
expensive  to  implement,  and  thus  methods  which  are  not  completely 
independent  of  the  director  system  must  be  considered  in  practice. 

A fourth  and  final  interpretation  feature  concerns  the 
relative  precision  of  the  target  state  vectors  obtained  from 
the  director  system  and  the  "true"  vectors  obtained  by  the  methods 
to  be  devised.  Thus  the  precision  required  of  the  "true"  vectors, 
in  order  that  they  may  be  considered  adequate  for  the  evaluation 
of  the  director  target  state  vectors,  is  not  specified. 

The  above  considerations,  which  may  lead  to  significantly 
different  interpretations  of  the  target  state  measurement  task, 
were  discussed  with  representatives  from  AFAL/RWT  and  from 
Honeywell,  and  some  decisions  were  made.  First,  it  was  decided 
that  the  director  system  target  state  vectors  would  be  evaluated 
at  a stage  close  to  the  final  output.  This  decision  allows  an 
assessment  of  the  effects  of  error  propagation  and  accumulation, 
although  information  on  the  effect  of  any  particular  stage  may 
be  lost.  In  a second  and  related  decision  it  was  determined 
that  the  target  state  vectors  would  be  obtained  in  the  attacker's 
body  frame,  since  this  is  the  coordinate  frame  in  which  these 
vectors  are  most  readily  available  at  states  close  to  the 
director  system  output. 
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A third  decision  determined  that  only  instrumentation 
planned  or  readily  available  for  tests  of  the  director  system 
would  be  used  to  devise  independent  methods  for  the  evaluation 
of  target  state  vectors.  Table  6 shows  the  real  and  test  inputs 
which  will  be  available  for  testing  the  director  system.  The 
real  inputs  could  be  used  in  an  actual  combat  encounter,  whereas 
the  test  inputs  are  available  for  test  purposes  only.  Note  that 
tracking  of  the  target  and  attacker  from  the  ground  or  from  an 
observation  aircraft  is  not  available  but  that  various  test 
inputs  transmitted  from  the  target  to  the  attacker  are  available. 
The  test  data  must  be  used  to  obtain  (for  the  evaluation  methods 
to  be  devised)  as  much  independence  from  the  director  system 
as  possible. 

In  a fourth  decision  it  was  determined  that  unproven  or 
unstable  analysis  techniques  would  be  avoided.  It  is  possible, 
for  example,  to  obtain  velocity  from  acceleration  if  the  initial 
integration  time  is  known,  and  it  may  be  possible  to  identify 
this  initial  time  as  the  time  at  which  integration  must  begin 
to  obtain  known  roll,  pitch,  and  yaw  angles  from  known  roll, 
pitch,  and  yaw  rates.  However,  this  is  an  apparently  unproven 
technique  which  may  require  an  unacceptable  amount  of  computer 
time  or  memory  if  it  were  implemented  on  the  attacker's  computer. 
Differentiation  to  obtain  velocity  from  position  or  acceleration 
from  velocity  is  an  example  of  a generally  unstable  technique  (in 
the  sense  that  initial  disturbances  are  accentuated) , although 
appropriate  smoothing  or  filtering  techniques  may  make  this 
procedure  acceptable. 
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TABLE  6 


REAL  AND  TEST  INPUTS 
TO  THE  DIRECTOR  SYSTEM 


REAL  INPUTS 


TEST  INPUTS 


ASCOT 

RADAR 

IRST 

HMS 

ISA 

ADC 

TIME  CODE 
MODE  CONTROL 


SSC 

MHC 

DME 

BARS 

TELEMETRY  FROM  TARGET 
ISA 
HARS 
ADC 


ASCOT 

- 

ADAPTIVE 

SCAN  OPTICAL  TRACKER 

IRST 

- 

INFRARED 

SEARCH  AND  TRACK 

HMS 

- 

HELMET  MOUNTED  SIGHT 

ISA 

- 

INERTIAL 

SENSOR  ASSEMBLY 

ADC 

- 

AIR  DATA 

COMPUTER 

DME 

- 

DISTANCE 

MEASURING  EQUIPMENT 

HARS 

- 

HEADING  ATTITUDE  REFERENCE  SYSTEM 

SSC 

- 

SOLID  STATE  CAMERA 

MHC 

- 

MINI  HUD 

CAMERA 

HUD 

- 

HEADS-UP 

DISPLAY 

T 


Finally,  it  was  decided  to  use  Monte  Carlo  methods  to 
obtain  error  probability  distributions  for  the  target  state 
vectors  calculated  according  to  the  interpretations  and  decisions 
discussed  above.  These  distribtuions  could  then  be  compared 
with  known  values  for  the  anticipated  accuracy  of  target  state 
vectors  from  the  director  system.  Anticipated  accuracy  is,  at 
present,  available  only  for  the  magnitude  of  the  acceleration  of 
the  target  relative  to  the  attacker:  a standard  deviation  of 

I 

about  Ig.  The  method  used  to  calculate  an  independent  acceleration  ■ 

‘ I 

and  an  associated  error  distribution  is  determined  by  the  i 

i 

decisions  and  interpretations  discussed  above,  and  it  involves 

the  use  of  simulated  ISA  and  HARS  data  from  both  target  and 

attacker.  The  mean  of  the  error  distribution  for  this  method, 

2 

in  the  most  realistic  simulation,  is  about  1.4  ft/sec  -much  less 
than  Ig.  Thus,  if  the  various  distributions  assumed  for  the 
Monte  Carlo  simulation  are  appropriately  conservative,  and  if 
the  errors  in  the  evaluation  method  are  sufficiently  uncorrelated 
with  the  director  system,  then  it  may  be  reasonable  to  conclude 
that  this  evaluation  method  is  adequate  for  the  verification  of 
target  acceleration  obtained  from  the  director  system.  Since 
anticipated  accuracy  figures  for  the  director  system  are  not 
available  for  velocity  and  position,  no  comparisons  can  be  made 
for  these  quantities. 

Note  that  no  attempt  is  being  made,  in  the  current  task, 
to  locate  the  source  of  any  discrepancies  which  may  be  observed 
between  the  target  state  vectors  obtained  from  the  director 
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system  and  the  "true"  vectors  obtained  by  independent  evaluation 
methods.  An  ability  to  identify  the  source  of  such  discrepancies 
may  form  an  important  extension  to  the  target  state  measurement 
task. 


3.2.2  Formulation  of  the  Monte  Carlo  Study 

The  objective  of  the  Monte  Carlo  study  is  to  obtain 
probability  distributions  for  the  errors  in  target  acceleration, 
velocity,  and  position  for  certain  evaluation  methods.  These 
methods  are  determined  by  the  decisions  made  concerning  the 
interpretation  of  the  task  objective,  as  discussed  above. 

Referring  again  to  Table  6,  none  of  the  real  inputs  should 
(ideally)  be  used  for  any  evaluation  method,  but  the  ISA  and  ADC 
inputs,  since  they  refer  only  to  the  attacker,  may  be  used  if 
necessary.  The  ISA  units  on  both  target  and  attacker  are  Honeywell 
H478  inertial  sensor  assemblies  and  are  strapped-down  acceler- 
ometer and  gyro  packages  whose  outputs  include  acceleration 
components  and  roll,  pitch,  and  yaw  rates.  True  air  speed  and 
angle  of  attack  are  available  from  the  Air  Data  Computer  on  each 
aircraft.  The  Heading  Attitude  Reference  System  is  part  of  the 
Stable  Coordinate  Reference  Group  (SCRGS)  which  in  turn  is  part 
of  the  MA-1  Fire  Control  System  on  both  target  and  attacker. 
(However,  HARS  is  not  a real  input  to  the  director  system) . The 
HARS  outputs  are  sines  and  cosines  of  roll,  pitch,  and  yaw  angles. 
The  Distance  Measuring  Equipment  consists  of  a transmitter  on 
the  attacker  and  a transponder  (remitter)  on  the  target.  The 
SSC  is  a Fairchild  CCD  solid  state  electronic  gunsight  camera 
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mounted  with  (but  independent  of)  ASCOT  on  the  attacker.  The 
KB-25  Mini  HUD  (Head-Up  Display)  camera  records  images  seen  by 
the  pilot  on  and  through  the  HUD,  including  images  of  the  target. 
Although  the  resolution  of  this  camera  may  be  high,  the  number 
of  frames  available  during  a simulated  combat  encounter  is 
probably  inadequate,  and  thus  a decision  was  mac.e  to  use  only 
the  SSC  to  simulate  line-of-sight  angle  data. 

Following  the  decisions  discussed  above  and  in  the  intro- 
duction, the  acceleration  of  the  target  relative  to  the  attacker 
is  available  only  using  ISA  and  HARS  data  from  both  target  and 
attacker.  Similarly,  the  velocity  of  the  target  relative  to 
the  attacker  is  avail-able  only  from  ADC  and  HARS  data  from  both 
target  and  attacker.  Finally,  the  position  of  the  target 
relative  to  the  attacker  is  available  only  from  DME  and  SSC 
data.  If  the  components  rather  than  the  magnitudes  of  the 
above  vectors  are  to  be  obtained,  then  the  appropriate  coordinate 
system  is  the  attacker's  body  frame. 

The  kinematics  of  the  above  methods  may  be  expressed  as 
follows.  Let  represent  the  vector  u of  the  object  x 

relative  to  the  object  y in  the  reference  frame  f;  e.g., 
is  the  velocity  of  the  target  relative  to  the  attacker  in  the 
attacker's  body  frame.  Let  T®  (q^)  represent  a coordinate 
transformation  from  the  initial  frame  s to  the  final  frame  e 
as  a function  of  the  quantities  q for  the  object  x;  e.g.,  T^ 

(<)>-»  0 , i);  ) is  a transformation  from  an  inertial  to  a body 

« a a 

frame  as  a function  of  the  roll((i)),  pitch  (6),  and  yaw  (ij;)  angles 
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of  the  attacker.  The  target  state  vectors  to  be  obtained  for 
comparison  with  the  same  vectors  from  the  director  system  are 


a a j a 

^t/a'  '^t/a'  ^ ^ P 


t/a’ 


</a  ” '’a'  ' 


a/i 


V 


t/a 


= M'^a'  «a'  ♦a>''t<'>t'  '*t'  »t> ' ’'v '“a*  Vi 


(21) 

(22) 


P?/a  = n'"a'  la’PU 


(23) 


Here  a,  v,  and  p are  acceleration,  velocity,  and  position  vectors; 
i,  V,  and  1 are  inertial,  velocity,  and  line-of-sight  frames; 
t and  a refer  to  target  and  attacker  or  to  the  corresponding 
body  frame;  <j)^,  0^,  and  are  roll,  pitch,  and  yaw  angles 
for  either  target  or  attacker  (for  x-a  or  t) , and  e and  n are 

3 3 

the  elevation  and  azimuth  of  the  target  relative  to  the  attacker. 
The  appropriate  transformation  matrices  are 


/ -se^  \ 

Cax  0 -Sa^\ 

0 1 0 I 

SUx  0 Ca^/ 
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(26) 


=l  C^xS\  ^^x  ^"x  S^x) 

\-Se  0 Ce„  / 

X X 

where  C and  S represent  sin  and  cos  and  ' indicates  matrix 
transposition.  Diagrams  of  the  methods  selected  for  the  inde- 
pendent evaluation  of  target  state  vectors  are  shown  in  Figures 


6,  7,  and  8. 

In  addition  to  the  kinem.atic  expressions,  probability 
distributions  for  the  occurrence  of  the  various  input  quantities 
are  required,  as  well  as  probability  distributions  for  the 
measurement  errors  in  these  quantities.  Tables  7 through  10 
show  these  input  quantities  and  their  estimated  occurrence  and 
accuracy  distributions.  Two  occurrence  distributions  for  the 
ISA  acceleration  components  (i.e.,  the  components  of  or 

were  tested.  The  first  is  a uniform  distribution  over  the 
range  -9g  to  +9g  for  each  component,  and  the  second  is  a much 
more  realistic  distribution  suggested  by  Capt.  Jerry  Kendrick 
of  AFAL/RWT:  the  x and  y components  of  acceleration  (in  the 
aircraft  body  frame)  are  normally  distributed  with  zero  mean 
and  standard  deviations  of  Ig  and  0.25g,  respectively,  and  the 
z component  is  distributed  according  to  the  log-  normal  distri- 
bution (8-4e®/2)g^  where  e~N(0,l).  This  distribution  has  a near 
cut-off  at  about  7.5g,  a mean  at  about  3.5g,  and  a "tail"  which 
extends  to  lower  accelerations.  Error  distributions  for  the 
ISA  acceleration  components  were  obtained  from  a table  of  instru- 
mentation signals  for  the  FI 06  provided  by  Capt.  Dave  Schoor 
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OCCURRENCE  AND  ACCURACY  DISTRIBUTIONS  FOR  THE 


OCCURRENCE  AND  ACCURACY  DISTRUBITONS  FOR 


OCCURRENCE  AND  ACCURACY  DISTRIBUTIONS  FOR 


r’ 


at  Tyndall  AFB.  Roll,  pitch,  and  yaw  ranges  and  measurement 
error  distributions  were  also  obtained  from  this  table.  Two 
different  occurrence  distributions  were  tested:  uniform 
distributions  over  the  permissible  ranges  and  a more  restricted 
piecewise  uniform  distribution  suggested  by  Ma j . Stewart  Cranston 
at  Tyndall  AFB.  Measurement  errors  for  computed  air  speed  and 
for  the  angle  of  attack  were  obtained  from  the  table  for  the 
F106  mentioned  above,  and  the  occurrence  distributions  for  these 
quantities  were  suggested  by  Capt.  Jerry  Kendrick  (the  side- 
slip angle  was  assumed  to  be  zero) . The  occurrence  and  error 
distributions  for  the  DME  were  suggested  by  Dr.  James  Yi  at 
Honeywell.  Finally,  the  required  occurrence  and  accuracy  distri- 
butions for  the  azimuth  and  elevation  angles  from  the  SSC  were 
estimated  from  information  provided  by  Capt.  H.  VJoodruff  and 
Mr.  Jerry  Watson  of  AFAL/RWT. 

The  procedure  for  the  Monte  Carlo  calculations  will  now 

t 3 V 

be  described.  Values  for  the  components  of  ^t/i' 

v^/i,  and  for  .he  angles  (j)^,  0^,  (p^,  6^,  a^,  e^,  and 

were  selected,  independently,  from  the  occurrence  distributions 

given  in  Tables  7 through  10.  Next  measurement  errors  for 

these  quantities  were  selected  from  the  accuracy  distributions 

a 


given  in  these  tables.  The  final  vectors,  ^t/a'  ^t/a' 


were  then  calculated  according  to  Eqs.  (21)  through  (26)  both 
with  and  without  the  addition  of  the  selected  measurement  errors 
to  the  selected  components  and  angles.  The  resulting  vectors 
were  subtracted,  and  the  entire  procedure  was  repeated  for  300 
selection  sets  for  each  vector  whose  error  distribution  was  to  be 
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determined.  Finally,  probability  distributions  of  the  differences 
between  the  square  root  of  the  sum  of  the  squares  of  the  com- 
ponents (for  vector  magnitude) , were  plotted.  The  final  result 
was  a set  of  simulated  error  distribu.>.  ns  for  the  evaluation 
methods  defined  by  Eqs.  (21)  through  (26)  and  diagrammed  in  Figures 
6,  7,  and  8.  These  evaluation  methods  are  proposed  as  independent 
methods  for  evaluating  the  accuracy  of  target  position,  velocity, 
and  acceleration  vectors  (relative  to  the  attaclcer  in  the  attacker's 
body  frame)  obtained  from  the  director  system. 

3.2.3  Results  of  the  Monte  Carlo  Stud 


Error  distributions  calculated  according  to  the  Monte 
Carlo  procedure  described  above  are  given  for  acceleration  in 
Figures  9 through  11.  These  distributions  refer  to  the  error 
in  the  acceleration  of  the  target  relative  to  the  attacker  in 
the  attacker's  body  frame  when  this  acceleration  is  obtained 
using  ISA  and  HARS  data  from  both  target  and  attacker.  Figure 
9 shows  the  magnitude  of  this  acceleration  error  calculated 
using  the  most  realistic  occurrence  distributions  for  the  ISA 
acceleration  components  and  HARS  angles  (log-normal  (a^,.) 

X/ 1 2 ' 

piecewise  uniform  (|)  and  9).  The  mean  of  this  distribution  is 

2 2 
about  1.4  ft/sec  - much  less  than  the  Ig  (32  ft/sec  ) anticipated 

accuracy  of  the  magnitude  of  the  acceleration  of  the  target 

relative  to  the  attacker  from  the  director  system.  In  addition. 

Figure  9 shows  that  the  probability  that  the  ISA-HARS  evaluation 

method  will  approach  an  error  of  Ig  is  very  small,  since  the 

"tail"  of  the  error  distribution  is  already  small  at  about  O.lg. 


Figure  10.  Error  Distribution  of  the  Magnitude  of  the  Acceleration 
of  the  Target  Relative  to  the  Attacker  Obtained  Using 
Simulated  ISA  and  HARS  Data 


Acceleration  Error,  ft/sec 


The  error  distributions  for  acceleration  magnitude  shown 

in  Figure  10  are  calculated  for  uniform  occurrence  distributions 

of  the  ISA  acceleration  components  (from  -9  to  +9g.)»  and  from 

uniform  or  piecewise  uniform  occurrence  distributions  for  (P , 9, 

and  ip.  These  error  distributions  have  means  considerably  larger 
2 

(about  3 ft/sec  ) than  the  distribution  of  Figure  9,  because  of 
the  much  greater  (and  probably  much  less  realistic)  probability 
of  large  acceleration.  However,  the  use  of  the  more  restricted 
piecewise  uniform  distribution  for  roll  and  pitch  angles  seems 
to  have  only  a minor  effect  on  the  final  error  distributions. 

Error  distributions  for  the  x,  y,  and  z components  of 
the  acceleration  of  the  target  relative  to  the  attacker  in  the 
attacker's  body  frame  are  shown  in  Figure  11.  These  distributions 
are  roughly  the  same  since  the  occurrence  and  error  distributions 
used  in  their  deviation  were  identical  (uniform  distributions 
from  -9  to  9g  for  the  ISA  acceleration  components  and  uniform 
distributions  for  6,  and  ip  over  their  allowed  ranges).  As  in 
Figures  9 and  10,  the  probability  that  any  error  in  acceleration 
will  exceed  Ig  is  very  small. 

Figures  12  and  13  show  error  distributions  for  the  magnitude 
of  the  velocity  and  position  vectors  of  the  target  relative  to 
the  attacker.  The  restricted  piecewise  uniform  occurrence 
distributions  for  ((),  9,  and  ip  were  used  for  the  velocity 
calculations.  The  mean  error  in  velocity  is  about  8 knots  and 
the  mean  error  in  position  is  about  4 ft.;  both  distributions 
are  skewed  toward  larger  magnitudes.  Although  these  distributions 
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Figure  12 . Error  Distribution  of  the  Magnitude  of  the  Velocity 

of  the  Target  Relative  to  the  Attacker  Obtained  Using 
Simulated  ADC  and  HARS  Data 


Velocity  Error,  Knots 


Figure  13.  Error  Distribution  of  the  Magnitude  of  the  Positi 
the  Target  Relative  to  the  Attacker  obtained  usin 
lated  DME  and  SSC  Data. 


Position  Error 
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provide  some  information  on  the  accuracy  of  the  ADC-HARS  method 
for  velocity  and  the  DME-SSC  method  for  position,  the  anticipated 
accuracy  of  the  director  system  for  these  vectors  is  not  known 
at  present,  and  thus  the  adequacy  of  these  methods  for  the 
evaluation  of  the  director  system  cannot  be  quantitatively  assessed. 

3.2.4  Suggestions  for  Future  Study 

As  discussed  in  the  introduction,  the  target  state  measure- 
ment task  has  been  interpreted  so  that  it  may  be  said  to  be  com- 
plete when  independent  methods  are  devised  and  shown  to  be  ade- 
quate for  the  evaluation  of  target  state  vectors  obtained  from 
the  director  system.  It  is  quite  likely  that  unacceptable 
discrepancies  will  be  found  between  the  target  state  vectors 
from  the  director  system  and  the  "true"  vectors  from  the  indepen- 
dent evaluation  methods.  In  the  current  task,  no  attempt  is 
made  to  locate  the  source  of  these  anticipated  discrepancies; 
however,  such  an  attempt  is  believed  to  be  essential  for  an 
adequate  evaluation  of  the  complete  director  system. 

There  are  at  least  three  apparent  approaches  to  the 
problem  of  "tracing  back"  observed  errors  in  the  director  target 
state  vectors.  First,  it  may  be  possible  to  perform  a sensitivity 

study  to  find  the  stages  in  the  director  system  at  which  major  | 

I 

discrepancies  originate.  Second,  the  target  state  measurement  ^ 

task,  as  interpreted  here,  might  be  carried  out  at  several  ' 

different  stages  in  the  director  system  (from  the  sensor  inputs  j 

to  the  HUD  output) . Finally,  a detailed  error  analysis  of  the  ■ 

j 

entire  director  system,  including  input,  processing,  and  output  i 
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errors,  could  be  made.  Results  from  any  of  these  approaches 
could  lead  to  significant  performance  improvements  in  the 
director  system. 


I 

; 3.3  IDENTIFY  AND  CORRECT  RADAR  LAG  PROBLEM 

3.3.1  Introduction 

I 3. 3. 1.1  Background 

I For  proposed  weapons  systems  to  be  effective, 

accurate  knowledge  of  the  target  state  is  needed.  The  EXPO  V and 

f 

1 

I Fire/Fly  programs  pointed  out  the  need  for  the  results  of  improved 

trackers  and  fire  control  algorithms.  To  this  end,  the  AFAL 

f initiated  the  Gunsight  Evaluation  Program  (Sight  Eval.)  to 

! 

evaluate  the  performance  of  existing  fire  control  systems.  This 
program,  conducted  at  Tyndall  AFB  using  an  instrumented  F-106 
which  tracked  towed  and  radio-controlled  targets  and  made  both 
live  and  dry  firing  missions,  produced  data  which  could  be  used 
for  the  analysis  of  existing  fire  control  systems. 

An  important  part  of  the  instrumentation 
aboard  the  F-106  was  the  MA-1  radar,  which  measured  range  and 
range  rate,  elevation  and  elevation  rate,  azimuth  and  azimuth 
rate.  It  is  the  purpose  of  this  report  to  analyze  the  performance 
of  the  t^A-l  radar  using  the  above  data. 

3. 3. 1.2  Motivation 

The  target  state  is  estimated  from  measurements 
of  range,  azimuth,  and  elevation,  along  with  the  first  and  second 
derivatives  of  range,  azimuth,  and  elevation.  The  derivatives 
are  obtained  from  suitable  algorithms  using  direct  measurements 
of  target  range  and  angular  position. 

The  MA-1  radar  provides  adequate  range 
information,  but  its  measurements  of  azimuth  and  elevation  are 
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' not  sufficiently  precise  to  be  useful  to  hand  off  target  ; 

1 ' 

position  to  an  electro-optical  sensor  (ASCOT),  whose  field  of 
view  is  much  narrower  than  the  MA-1,  and  which  measures  the 
angular  position  of  the  target  with  sufficient  accuracy. 

The  major  system  integration  problem  with 
the  above  technique  is  that  the  two  fields  of  view  of  the  ASCOT 
are  limited  to  + 17  mrad  and  +43  mrad,  whereas  the  MA-1  reports 
angular  positions  which  are  in  error  by  as  much  as  + 50  mrad.  j 

It  is  the  purpose  of  this  analysis  to  improve  the  radar 
estimates  of  the  angular  position  such  they  have  a high  proba- 
bility of  lying  within  the  + 17  mrad  window  of  the  ASCOT. 

3. 3. 1.3  Scope  and  Objectives  of  the  Problem 

There  are  two  different  methods  which  can  be 
used  to  do  the  above.  The  first,  called  the  hardware  approach, 
is  to  modify  the  design,  mounting,  or  calibration  of  the  MA-1 
such  that  the  desired  degree  of  accuracy  can  be  achieved.  The 
second,  a software  approach,  is  to  look  for  a suitable  algorithm 
which  will  produce  correction  terms  such  that  the  revised  radar  | 

estimates  lie  within  the  allowable  error  tolerances. 

The  principal  drawback  to  the  hardware  approach 
is  the  lack  of  access  to  the  design  and  installation  data  on  the 
modified  MA-1  system  and  to  the  schedules  of  the  alignment  and 

I 

calibration  procedures  performed  on  the  system  during  the  Sight  ' 

Eval  tests.  j 
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The  second  approach  is  a possible  alternative. 


The  ROLM  16/64,  a rugged,  high-speed,  airborne  computer  is  a 
component  of  the  proposed  fire  control  system,  and  could  con- 
ceivably contain  provisions  for  correcting  the  angular  data 
provided  by  the  MA-1.  It  is  along  these  lines  that  this  study 
has  been  directed. 

Two  distinct  efforts  related  to  a software 
approach  were  pursued.  One  effort  involved  a relatively  straight- 
forward graphical,  curve-fit  attack  attempting  to  relate  tracking 
error  to  tracking  rate  functions.  The  second  effort  involved 
an  analysis  of  the  relative  angular  rate  of  change  of  the  Line 
of  Sight  (LOS)  between  the  aircraft  body  axis  and  the  target 
followed  by  a similar  analysis  of  angular  rates  between  aircraft 
and  target  relative  to  an  inertial  coordinate  system. 

3.3.2  Data  Collection 

3. 3. 2.1  General 

The  original  purpose  of  the  Sight  Eval  pro- 
gram was  to  test  and  evaluate  fire  control  systems,  especially 
the  HLGS  (Hot  Line  Gun  Sight)  and  the  DALCOS  (Digital  Advanced 
Lead  Computing  Optical  Sight) . As  mentioned  previously,  the 
objective  of  this  study  is  to  determine  the  magnitude  and  causes 
of  the  angular  position  errors  in  the  angular  position  measurements 
of  the  MA-1  Fire  Control  System,  and  to  develop  an  algorithm  to 
reduce  these  errors.  Of  primary  concern  is  the  difference 
between  the  "true"  target  position,  determined  from  a 16  mm 
film  record,  and  the  position  reported  by  the  MA-1  radar. 
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The  radar  target  position  and  rates  are  determined  by  resolvers 
mounted  on  the  MA-1  antenna  gimbals,  and  by  interceptor  rate 
gyros  and  radar  range  measurements. 

The  data  available  for  this  study  consisted 
of  analog  signals  recorded  on  magnetic  tape  and  target  angular 
position  recorded  on  16  mm  film.  After  each  mission,  the  magnetic 
tape  and  film  data  were  converted  to  digital  form  and  recorded 
on  tape. 

3. 3. 2. 2 Description  of  Test  Missions 

Ninety-five  live  firing  missions  were  planned 
and  executed  for  the  purpose  of  evaluating  gunsight  algorithms, 
but  only  39  of  these  missions  were  recorded. 

Figure  14  shows  the  relationship  between 
target  and  attacker  in  a typical  firing  encounter.  In  one 
mission,  the  target  was  remotely  piloted  PQM-102.  In  the 
remaining  missions,  the  targets  were  towed.  The  towed  targets 
were  either  the  30  ft  long  fiberglass  FIGAT  or  the  15  ft  long 
aluminum  DART. 

The  firing  tests  were  conducted  at  different 
ranges,  altitudes,  angle-off  at  firing,  and  target  acceleration. 
The  target  was  in  either  a Ig,  2g,  or  4g  turn.  The  attacker 
coordinated  his  pursuit  so  that  firing  began  at  prescribed  angles- 
off  (10°,  20°,  30°,  or  more)  and  ranges  (1000,  1500,  2000,  or 
3000  ft)  at  altitudes  of  10,000,  15,000,  and  20,000  ft.  Figure 
15  is  a matrix  showing  the  combinations  of  these  variables  which 
were  used  in  the  tests. 
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Figure  14  Firing  Profile  for  Actu^ 
VJith  a Drone  Target  (PQ! 


3. 3. 2. 3 Instrumentation  and  Recorded  Data 


A block  diagram  showing  the  transfer  of  analog 
data  from  the  sensors  to  the  recording  tape  is  shown  in  Figure 
16.  The  ISA  (H  478  Inertial  Sensor  Assembly)  measures  ownship 
acceleration  in  rectangular  components.  The  IRP  (Inertial 
Reference  Package)  gives  ownship  pitch,  yaw,  and  roll  rates. 

The  SCRG  (Stable  Coordinate  Reference  Group)  - AHARS  (Altitude 
Heading  and  Reference  System)  gives  ownship  pitch,  yaw,  and  roll. 

Radar,  operating  in  conjunction  with  the  CADC 
(Central  Air  Data  Computer) , gives  range  and  range  rate.  The 
IR,  and  infra-red  tracker,  provides  passive  angular  target 
direction  for  powered  targets.  These  signals  are  input  to  an 
HCM-204  Fire  Control  Computer.  The  3 -vane  is  a yaw  attitude 
sensor. 

All  the  signals  described  above  are  fed  into 
a signal  conditioner,  whose  function  is  to  adjust  the  amplitude 
of  the  incoming  analog  signals  to  a 0-5  volt  level,  and  to 
generate  a Pulse  Amplitude  Modulated  (PAM)  train,  which  is  an 
ordered  sequence  of  all  signals.  The  sets  of  data  on  the  PAM 
train  are  separated  by  about  70  msec.  The  output  of  the  signal 
conditioner  is  recorded  on  a Genesco  tape  recorder.  The  output 
of  the  IRIG  - B Time  Code  Generator  is  recorded  directly  onto 
the  tape,  to  the  nearest  thousandth  of  a second. 

The  conversion  of  data  from  the  ISA  and  the 
MA-1  sensors  into  a PAM  train  by  the  signal  conditioner  is 
shown  in  Figure  17.  The  PAM  train,  together  with  ballistic 
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information  for  the  type  of  projectile  used,  are  fed  into  a HDC- 
301  computer.  Using  an  algorithm  peculiar  to  each  of  the  gunsights 
to  be  compared,  the  computer  determines  the  future  azimuth  and 
elevation  of  the  target.  An  interface  box  converts  these  signals 
into  a form  such  that  they  may  be  rendered  as  symbols  on  the 
HUD  (Heads  Up  Display) . They  are  simultaneously  seen  by  the 
pilot  and  recorded  on  flight  film  by  a KB  25  or  KB  26A  camera. 

The  conversion  of  the  flight  film  and  flight 
tape  to  digital  form  and  the  result  recorded  on  tape  in  a form 
compatible  with  the  CDC  6600  is  shown  in  Figure  18.  The  empty 
box  symbolizes  unknown  processing  and  transcribing  operations 
which  occurred  in  some  cases.  The  result,  the  WPAFB  tapes, 
were  the  raw  material  for  this  analysis. 

3.3.3  Preliminary  Data  Evaluation  and  Processing 
3. 3. 3.1  Description  and  Selection  of  Data 

The  information  on  the  WPAFB  tapes  came 
from  data  recorded  during  test  missions  flown  at  Tyndall  AFB 
during  1975  and  1976.  Figure  19  lists  the  missions  for  which 
data  were  recorded  and  supplied  to  WPAFB.  Each  mission  consisted 
of  a number  of  passes,  and  only  data  from  wet  (firing)  runs  were 
recorded.  A study  of  the  data  led  to  the  classification  of  the 
missions  into  seven  groups,  two  of  which  contain  but  one  mission. 
Data  from  each  pass  were  considered  in  the  light  of  pilot  comments. 
In  some  cases  of  obvious  malfunction;  e.g.,  no  radar  lock-on, 
the  pass  was  disregarded. 
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Recorded  Missions 
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Each  pass  record  consists  of  the  PAM  data 
plus  signals  from  the  time  signal  generator.  Thirty-seven 
channels  of  data  recorded  on  the  WPABF  tapes  are  listed  in 
Table  11. 

Data  relevant  to  this  analysis  were  the  time, 
range,  pitch  rate,  yaw  rate,  radar  antenna  azimuth,  radar  antenna 
elevation,  film  target  position  azimuth,  and  film  target  position 
elevation.  When  the  validity  of  the  data  was  in  doubt,  acceler- 
ometer, rate  gyro,  and  heading  information  were  examined  as 
possible  sources  of  error. 

3. 3. 3. 2 Evaluation  and  Preliminary  Processing  of  Data 

The  film  records  were  made  on  16  mm  film  by 
a KB-25  or  KB-26A  camera  operating  at  24  frames  per  second. 

Air  Force  engineers  have  estimated  a maximum  error  of  1 mrad 
due  to  film  sensitivity.  Alignment  has  reduced  parallax  to  about 
the  same  magnitude.  Because  of  the  finite  size  of  the  target 
as  seen  on  the  film,  a point  on  the  target  - usually  the  tail  - 
was  chosen  as  the  target  location.  However,  inconsistency  in 
film  reading  practice  could  have  led  to  variations  of  as  much 
as  5 mrad  in  the  target  position. 

At  close  range,  an  error  of  from  10  to  20 
mrad  could  have  been  introduced  due  to  differences  between 
the  point  on  the  target  selected  from  film  reading  and  that 
which  was  most  strongly  reflecting  the  radar  signal. 

Film  elevation  and  azimuth  data  and  radar 
range  data,  elevation,  and  azimuth  data  were  smoothed  by  means 
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TABLE  11 


TAPE  DATA  VARIABLES 


WORD  NO.  IN  RECORD 


VARIABLE 


UNITS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40-50 


IRIG  TIME 

TRACKER  MODE(0=Range,  1=IR) 
MACH  NUMBER 
IMPACT  TEMPERATURE 
STATIC  PRESSURE 
ANGLE  OR  ATTACK 
RANGE 

RANGE  RATE 
MANUAL  RANGE 
TRIGGER  #1  (l.=  ON) 
TRIGGER  #2  (l.=  ON) 

RANGE  ON  TARGET  (l,=LOCK) 
DELAYED  RANGE  ON  TET  (l.=LOCK) 
ROLL  RATE 
PITCH  RATE 
YAW  RATE 
X ACCELEROMETER 
Y ACCELEROMETER 
Z ACCELEROMETER 
SIDESLIP  ANGLE 
EVENT  (l.=  ON) 

HEADING 

ROLL 

PITCH 

RADAR  ANTENNA  AZIMUTH 
RADAR  ANTENNA  ELEVATION 


MSEC 


2 

LES/FT 

RAD 

FT 

FT/SEC 

FT 


RAD/SEC 
RAD/SEC 
RAD/SEC „ 
FT/SECp 
FT/SEC2 
FT/SEC 
RAD 


RAD 

RAD 

RAD 

RAD 

RAD 


IR  GIMBAL,  AZIMUTH  RAD 

IR  GIMBAL,  ELEVATION  RAD 

AMMO  TEMPERATURE  (UNCALIB.) 

FILM  TARGET  POSITION,  AZIMUTH  MRAD 

FILM  TARGET  POSITION,  ELEVATION  MRAD 

FILM  LCOS  POS,  AZIMUTH  MRAD 

FILM  LCOS  POS,  ELEVATION  MRAD 

FILM  HLES  POS,  AZIMUTH  MRAD 

FILM  HLES  POS,  ELEVATION  MRAD 

FILM  DESIGNATOR  POS,  AZIMUTH  MRAD 

FILM  DESIGNATOR  POS,  ELEVATION  MRAD 

(BLANK)  
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of  a polynomial  function,  although  the  raw  film  data  were 
reasonably  smooth. 


The  rms  error  between  the  radar  range  data 
and  the  fitting  function  was  calculated  for  90  passes  and  the 
results  plotted  in  Figure  20.  The  rms  variation  in  radar  range 
lies  between  12  and  53  ft  (4  to  16  m)  and  shows  no  strong  de- 
pendence on  range.  If  one  assumes  that  the  maximum  rms  deviation 
of  53  ft  corresponds  to  3 a , then  the  noise  in  the  radar  range 
measurements  has  a standard  deviation  of  less  than  6 m. 
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Figure  20  RMS  Range  Error 
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3.3.4  Analysis  and  Correction  of  Lag  Problem 
3. 3. 4.1  Linear  Regression 

Preliminary  plots  of  radar  and  film  data 
indicated  that  the  radar  data  was  usually  lagging  the  film  data. 
To  use  the  radar  to  predict  the  target  position  at  a future 
time  point  or  to  aim  the  ASCOT,  it  is  necessary  to  compensate 
for  the  lag.  The  ideal  system  adjustment  would  be  a software 
correction  if  a functional  relationship  could  be  identified 
which  consistently  reduced  the  radar  lag.  This  functional 
relationship  also  had  to  reduce  the  radar  lag  to  + 15  MRA.D,  the 
ASCOT  lock-on  constant.  Average  radar  lag  was  calculated  by 
subtracting  film  position  from  radar  position  for  each  pass. 

It  was  noted  that  the  average  radar  lag  varies  considerably 
between  passes  on  the  same  mission. 

If  it  is  assumed  that  radar  tracking  error 
is  a function  of  line-of-sight  rate,  the  general  form  of  the 
error  is: 

= Af  + + C^q 

=^2*  ®2"r  + <^2'^ 
where  E - Error 

e - Target  elevation  rate  of  change  with  respect  to 
interceptor  airframe  axis 

Target  azimuth  rate  of  change  with  respect  to 
interceptor  airframe  axis 

q - Interceptor  pitch  angle  rate 

r - Interceptor  yaw  angle  rate 

A,B,C,  - Constants  derived  by  regression  program 

The  target  crossing  and  consequent  radar 
tracking  error  was  assumed  to  be  a combination  of  a bias.  A, 
plus  a constant,  B,  times  the  azimuth  or  elevation  tracking  rate 
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Tracking  Rates  (!!„»£„)  vs  Mean  Radar  Error  (y) 
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of  the  radar  antenna,  plus  a constant,  C,  times  the  pitch  or 
yaw  angle  rates  of  the  interceptor.  The  missions  were  divided 
into  seven  chronological  groupings  and  a regression  analysis 
was  performed  on  each  group  to  determine  the  coefficients  A,  B, 
and  C.  The  data  groupings  (as  defined  in  Section  3. 3. 3.1)  are 
shown  in  Figures  21  thru  27. 

The  results  of  the  regression  indicated  that 
the  significant  contribution  was  always  the  term.  This  is 

not  surprising  as  the  pitch  and  yaw  rates  were  small  compared 
with  e as  was  the  C coefficient.  The  total  contribution  of  Cq 
or  Cr  was  thus  negligible.  The  A coefficient  defines  the  bias 
in  the  radar.  The  bias  contribution  while  small,  varied  between 
data  groups  indicating  a negligible  constant  bias.  Ideally  if 
the  radar  lag  is  a true  function  of  crossing  rate,  a zero  crossing 
rate  implies  no  radar  lag.  For  these  reasons  the  functional 
relationship  for  all  missions  was  simply  expressed  as 


B, 

; B = 

-135 

e 

1 R 

1 

E = 

: B^  = 

-250 

n 

2 R 

2 

The  mean  radar  error  components  were  then 
calculated  for  each  pass.  The  tracking  error  correction  was 
applied  and  the  mean  radar  error  was  again  calculated.  Histograms 
of  the  mean  radar  errors  are  shown  in  Figures  28  thru  31.  The 
mean  error  plotted  can  be  expressed  as 

- e - E = E 
R F e ec 


■n  - ri  - E 
R 'f  q 


where 


= E 


Tic 

Radar  Elevation 
Film  Elevation 


Tlj^  Radar  Azimuth 
Eg  Elevation  Tracking  Corrector 
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Figure  21  Data  Group  I 
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Figure  27  Data  Group  VII 


Figure  28  Average  EL  Error  Before  Tracking  Correction 


Figure  29  Average  EL  Error  After  Tracking  Correction 


Figure  31  Average  AZ  Error  After  Tracking  Correction 


Azimuth  tracking  Correction 


1 


E„  Elevation  error 
Gc 

E„  Azimuth  error 

' Ic 

The  tracking  error  correction  did  reduce  the 
radar  mean  lag  for  a number  of  the  passes.  However,  a significant 
number  of  the  passes  were  still  greater  than  the  17  MRAD  con- 
straint. This  suggests  that  error  sources  other  than  crossing 
rates  are  contributing  to  the  radar  lag. 
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3. 3. 4. 2 Line-of-Sight  Rate  of  Change  Analysis 

Validity  of  Sight  Eval  Data 

To  provide  an  alternative  to  graphical  methods  of  studying 
the  radar  lag  problem,  a program,  SDATA,  was  written  to  reduce  the 
Sight  Eval  data.  SDATA  first  smooths  the  Sight  Eval  data  and 
then  uses  an  Eulor  parameter  integration  scheme  to  compute  the 
principal  variables  during  each  engagement.  The  important 
variables  computed  are  range,  range  rate,  pursuer  and  target 
positions  and  velocities,  the  absolute  LOS  turn  rate  and  the  LOS 
turn  rate  relative  to  the  pursuer,  pursuer  turn  rate,  headings 
of  both  pursuer  and  target. 

The  program  SDATA  had  a second  purpose  which  was  to 
determine  if  the  Sight  Eval  data  was  invalid.  Once  the  SDATA 
integration  scheme  is  initialized,  it  is  possible  to  compare  the 
computed  values  of  pursuer  roll,  pitch,  and  yaw  angles  and  pur- 
suer speed  with  those  values  obtained  from  smoothed  Sight  Eval 
data.  These  comparisons  indicate  that  the  Sight  Eval  data, 
although  much  noisier  than  would  be  hoped,  does  furnish  a reason- 
able basis  for  a study  of  the  radar  lag  problem,  assuming,  of 
course,  that  none  of  the  scaling  or  other  operations,  performed 
on  the  tapes  prior  to  our  effort,  introduced  significant  error. 
Documentation  has  not  been  obtained  on  some  aspects  of  these 
operations.  Two  other  qualifications  are  in  order.  The  geometry 
of  most  of  the  engagements  or  passes  in  the  Sight  Eval  data  are 
very  similar  in  geometry.  This  raises  the  question  of  the  con- 
fidence one  could  have  in  a derived  fit  for  a fundamentally 
different  engagement  geometry.  Another  problem  is  that  most 
passes  are  of  short  duration,  often  less  than  a second,  rarely 
over  2 seconds.  Would  a fit  be  valid  for  larger  periods  of  time? 

It  was  hoped  that  an  analysis  of  a large  number  of 
engagements  would  allow  the  derivation  of  an  equation  giving 
radar  lag  angle  as  some  function  of  the  various  engagement  turn 
rates.  Three  turn  rates  that  seem  likely  candidates  are  true 


LOS  turn  rate  with  respect  to  inertial  coordinates,  pursuer 
ownship  turn  rate,  and  the  relative  LOS  turn  rate  with  respect  to 
the  pursuer. 


i 


With  this  goal  in  mind,  a number  of  fitting  equations 
were  derived  and  added  to  SDATA.  A large  number  of  engagements 
were  then  analyzed  in  SDATA  and  the  results  of  each  of  the  fits 
for  each  engagement  was  tabulated.  Two  different  types  of  fit 
were  utilized  and  their  deviation  and  intent  will  now  be  outlined. 

In  the  following  discussion  all  turn  rates  and  the  vector 
radar  lag  angle  itself  are  expressed  in  terms  of  their  components 
in  the  SDATA  inertial  x y z coordinate  system. 

Radar  lag  angle  is  represented  as  a vector  as  follows: 


I 


i 

I 

♦ 


In  Figure  32  rn  and  rl  are  unit  vectors  along  the  radar  axis 

and  along  the  true  line  of  sight  to  the  target  respectively,  both 

vectors  emanating  from  the  pursuing  aircraft.  Then  by  definition 
— m X rl 

9 = 0 1 21: ^l-  This  procedure  is  valid  since  9 is  of  relatively 

' rn  X rl ' 

small  magnitude  and  it  has  the  advantages  of  eliminating  many 
annoying  sign  convention  problems  as  well  as  simplifying  the 
formulation  of  fit  equations. 

Let  wa  be  the  true  LOS  turn  rate,  wb  the  pursuer  ownship 
turn  rate,  and  ^ the  angular  rate  of  change  of  the  LOS  with 
respect  to  the  pursuing  aircraft.  Let  the  number  of  time  points 
in  an  engagement  be  n.  The  first  fit  equations  are  of  the  type 

a.  = a + b + c 

1111 

where  the  constants  a,  b,  and  c are  chosen  so  as  to  minimize 


a . = c wc . 

1 1 

The  second  type  of  fit  tried  included  a bias  vector  term  b which 
was  intended  to  eliminate  variations  between  engagements  caused 
by  differing  initial  conditions  at  radar  lockon.  This  equation 
is  of  the  form 

^ = a wa . + b 
1 1 

where  both  the  constant  a and  the  three  constant  components  of  b, 
bx,  by,  and  bz  must  be  found.  Again,  setting 


M ^ 6E  ^ ^ ' 

Ta  6bx  ^y 

allows  a straightforward  computation  of  these  constants. 

To  obtain  a measure  of  the  reduction  in  lag  angle  for 
a given  fit,  a sample  mean  and  a sample  standard  deviation  : 

calculation  is  performed  on  both  the  lag  angle  6^  and  on  the  i 

"unexplained"  part  of  lag  angle  follows:  , 


fit  = - Z |0.-a. 
av  n ' 1 1 


_ _ 2 j 

Z I 0 . -a . I - n f it 
i=l  ' ^ ^ 


Fit-Equation  Validation 


Each  of  the  above  fits  reduces  the  average  lag  angle, 
sometimes  quite  dramatically,  but  of  course  the  key  question  is 
"Is  there  a fit  whose  constants  remain  relatively  stable  from 
engagement  to  engagement?"  Apparently,  there  is  not.  Each  of 
the  fits  tried  exhibits  a randomness  in  its  constants  over  a 
number  of  passes  even  within  the  same  mission. 

Of  the  fits  examined,  the  two  giving  the  best  overall 
reduction  in  the  mean  radar  lag  angle  while  at  the  same  time 
having  the  least  variation  in  the  fitting  constants  were 
cr=a  wa^,  the  fit  to  lag  angle  as  a constant  times  the  true  LOS 
turn  rate,  and  cr=c  wc^,  the  fit  as  a constant  times  the  LOS 
turn  rate  relative  to  the  pursuing  aircraft.  The  following  data 
are  based  on  106  passes  from  the  latest  tape  data.  No  pass  was 
included  that  had  a mean  lag  angle  in  excess  of  90  mils.  The 

average  lag  angle  over  these  passes  was  34.8  mils  with  an  average 

standard  deviation  of  7.45  mils. 

Reference  Table  13,  list  of  passes  checked,  using  the 
fit  a^=a  wa^  gave  an  average  fitted  lag  angle  of  11.8  mils  with 
an  average  standard  deviation  of  5.04  mils.  But  the  constant 
a varied  between  0.028  to  0.432  with  the  bulk  of  the  a's  falling 
between  0.10  and  0.26.  The  c's  varied  between  the  extremes  of 
0.140  and  0.325.  This  large  variation  in  fit  constants  shows 
clearly  the  difficulty  of  attempting  to  correct  for  lag  angle 
in  this  manner.  This  variation  in  the  two  best  performing  fits 
tried  is  a result  of  the  fact  that  large  lag  angles  occur  with 

both  high  and  low  turn  rates  as  do  the  small  lag  angles.  Assuming 
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TABLE  13  LIST  OF  PASSES  CHECKED 


Fit 

av 

Fit 

av 

fit  . 
sd 

Mission 

Pass 

3( sec) 

c(  sec) 

a 

av(mils) 

3(mils) 

C(mils ) 

0 

sd{mils) 

s(mils) 

7(mils) 

76077 

1 

0.072 

0.146 

10.7 

4.0 

5 . 6 

2.3 

2.2 

2.9 

" 

2 

0.108 

0.187 

14.9 

5.2 

5 . 6 

2.3 

2.0 

2.5 

!• 

3 

0.125 

0.212 

21.6 

6.1 

6.0 

3.3 

3.3 

1.7 

*• 

4 

0.029 

0.044 

3.3 

5 . 5 

5.8 

4.0 

2.8 

3.1 

*' 

3 

0.101 

0.177 

30.2 

14.9 

7.  3 

1.6 

3.4 

2.2 

" 

6 

0.048 

0.068 

15.9 

3.1 

4.2 

3.0 

2.2 

3.6 

76073 

1 

0.051 

0.129 

8.3 

5.7 

5.6 

1.8 

2.1 

2.2 

2 

0.114 

0.220 

23.4 

7.3 

5.7 

3.8 

3.2 

3.5 

3 

0.041 

0.294 

15.9 

14.9 

11.3 

4.5 

5.1 

4.7 

4 

0.138 

0.419 

17.5 

3.4 

6.4 

3.1 

1.5 

4.5 

" 

5 

0.029 

0.215 

18.8 

18.3 

15.3 

13.7 

13.9 

11.3 

6 

0.159 

0.228 

47.5 

9.3 

8.0 

8.9 

4.7 

2.6 

" 

7 

0.120 

0.180 

35.9 

17.5 

12.5 

3.2 

5.0 

3.1 

" 

8 

0.170 

0.223 

45.0 

15.3 

8.6 

1.9 

4.0 

3.9 

76093 

1 

0.130 

0.183 

37.5 

•<  . X 

9.2 

3.7 

1.4 

1 . ■’ 

*» 

2 

0.113 

0.140 

30.2 

15.8 

14.2 

4.8 

2.8 

2.0 

II 

3 

0.120 

0.150 

33.7 

1.7 

4.9 

2.0 

1.2 

0.5 

M 

4 

0.158 

0.217 

37.6 

15.2 

13.7 

2.3 

1.0 

1. 1 

76085 

1 

0.054 

0.555 

26.5 

24 . 4 

13.6 

4.7 

3.  5 

6 . 1 

76092 

i_ 

0.240 

0.243 

35.0 

21.1 

29 . 5 

14.7 

13.7 

12.1 

76093 

1 

0.082 

0.365 

25.1 

22.0 

11.3 

7.4 

9.0 

4.9 

76110 

1 

0.432 

0.728 

87.1 

50.5 

59.2 

27.2 

21.7 

19.3 

If 

4 

0.366 

0.561 

82.2 

41.0 

29.6 

39.0 

14.4 

12.5 

II 

6 

0.422 

0.851 

78.9 

41.8 

19.5 

39.4 

18.6 

9.3 

76111 

1 

0.134 

0.161 

23.2 

5 . 5 

5.3 

4.0 

2.2 

2.8 

" 

2 

0.171 

0.246 

21.8 

7.1 

3.2 

1.6 

1.5 

1.3 

3 

0.145 

0.202 

40.6 

12.9 

4.6 

3.6 

2.4 

1.4 

rl 

4 

0.130 

0.180 

36.8 

6.8 

4.5 

5.9 

3.4 

2 . 1 

" 

5 

0.149 

0.229 

45.8 

14.7 

2.0 

5.0 

1.3 

o'.  6 

76113 

1 

0.188 

0.243 

28.8 

3.9 

4.2 

2.6 

2.1 

1.  8 

" 

2 

0.170 

0.191 

30.3 

6.5 

9.8 

1.5 

3.4 

3.3 

•* 

3 

0.184 

0.258 

42.0 

8.5 

5.6 

4,5 

3.6 

3.0 

If 

4 

0.067 

0.034 

15.4 

5.4 

5.2 

2.5 

3.4 

3.5 

5 

0,174 

0.207 

33.6 

3.8 

7.2 

1.8 

1.6 

2.0 

6 

0.251 

0.301 

74.9 

14.7 

8.7 

12.1 

4.7 

4.7 

•• 

7 

0.221 

0.251 

66.5 

12.1 

14.3 

2.9 

6.7 

1.2 

76119 

2 

9.120 

0.306 

11.8 

4.7 

4.4 

3.  3 

2.0 

3.1 

" 

3 

0.114 

0.168 

16.5 

5.1 

7.2 

3.7 

2.8 

4.7 

4 

0.166 

0.278 

25.8 

4.3 

7.5 

1.9 

1.3 

3.2 

3 

0.097 

0.312 

20.6 

16.3 

14.0 

6.5 

8.0 

5.6 

6 

0.201 

0.291 

36.4 

7.6 

17.0 

2.7 

3.2 

10.2 

7 

0.223 

0.306 

45.9 

9.7 

7.2 

3.0 

6.3 

2.7 

76121 

1 

0.037 

0.088 

5.7 

4.9 

4.7 

3.2 

2.4 

2.3 

If 

2 

0.069 

0.160 

11.8 

9.8 

9.7 

3.0 

2.5 

1.5 

3 

0.080 

0.243 

10.7 

7.3 

6.4 

1.6 

1.  5 

2.1 

If 

4 

0.079 

0.215 

14.6 

7.2 

6.7 

4.5 

4.7 

4.8 

It 

3 

0.236 

0.282 

50.3 

4.5 

9.8 

3.5 

1.  2 

1.1 

6 

0.132 

0.163 

32.3 

8.3 

3.9 

2.4 

1.9 

0.8 

7 

0.119 

0.198 

24.1 

6.1 

13.9 

2.6 

1.6 

2.0 

If 

3 

0.072 

0.075 

17.6 

12.4 

7.5 

1.9 

2.8 

2.8 

9 

0.202 

0.256 

55.1 

7.1 

5.8 

5.7 

2.5 

3.0 

76126 

1 

0.066 

0.210 

12.2 

3.3 

13.7 

1.5 

1.9 

1.4 

2 

0.135 

0.179 

24.3 

18.0 

7.4 

1.5 

3.4 

1.8 

** 

3 

0.113 

0.246 

11.5 

7.6 

13.3 

3.9 

2.3 

2.9 
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TABLE 

13 

LIST  OF 

PASSES 

CHECKED 

(CONTINUED) 

Fit 

fit.,. 

Fit  , 

fit  , 

av 

av 

sd 

sd 

Mission 

Pass 

a(sec) 

c(sec) 

0 

av(mils)  aimils) 

c(mils) 

9 

sdfnils) 

avmilsj 

c(inils) 

76131 

1 

0.115 

0.243 

19.4 

18.2 

13.3 

21.4 

15.2 

15.  7 

" 

2 

0.042 

0.050 

14.9 

14.2 

14.5 

2.1 

2.6 

2.5 

" 

3 

0.166 

0.302 

47.  3 

14.1 

18.2 

17.1 

6 . 2 

6.9 

4 

0.  105 

0.253 

24.1 

14.7 

11.6 

5.0 

3.1 

3.3 

" 

5 

0.107 

0.171 

28.4 

9.3 

6.0 

2.9 

2.2 

2.1 

tl 

6 

0.059 

0.206 

12.  1 

3.3 

6.4 

2.8 

2.8 

2.9 

" 

7 

0.194 

0.305 

58.3 

19.5 

6.3 

5.8 

1.4 

4.6 

" 

3 

0.229 

0.345 

72.5 

14.2 

18.0 

12.1 

4.4 

2 . 6 

10 

0.059 

0.091 

21.6 

15.4 

15.8 

4.  3 

2.3 

2.0 

76162 

1 

0.065 

0.243 

31.0 

76.7 

71.0 

60.0 

63.5 

65.5 

" 

2 

0.143 

0.273 

44.1 

22.7 

15.1 

16.7 

4.0 

2.2 

3 

0.151 

0.205 

40.0 

14.7 

12.6 

5.2 

3.5 

2.  3 

4 

0.120 

0.231 

36.8 

21.3 

20.2 

23.  9 

27.6 

13.0 

5 

0.198 

0.243 

56.3 

16.9 

18.1 

3.7 

1.9 

1.4 

5 

0.151 

0.326 

61.2 

32.5 

14.1 

15.3 

1.  3 

1.7 

7 

0.194 

0.225 

65.4 

54.5 

49.3 

53.7 

31.2 

63.3 

76163 

1 

0 . 250 

0.386 

34.5 

12.4 

11.3 

2.3 

4.7 

4.3 

'* 

3 

0.119 

0.209 

36.1 

9.7 

3.1 

2.3 

3.  3 

1 . 5 

'* 

4 

0.249 

0.306 

49.0 

9.8 

10.1 

3.2 

3.7 

3.1 

tt 

5 

0.293 

0.376 

42.3 

9.7 

10.2 

4.2 

4.0 

2.9 

" 

6 

9.169 

0.201 

51.3 

9.2 

5 . 1 

4.3 

“ 2 

3.1 

'* 

7 

0.069 

0.  146 

14  . 5 

11.1 

10.0 

1 

2 , 3 

■>  _ 2 

" 

8 

0.251 

0.312 

52.3 

8.2 

4 . 3 

4.3 

4.7 

3.4 

'* 

9 

0.260 

0.235 

43.9 

10.4 

13.2 

5.4 

3.8 

3 . 1 

76166 

1 

0.031 

0.174 

13.8 

13.5 

10.0 

6.2 

5.6 

4.4 

" 

2 

0.032 

0.183 

5.7 

4.9 

3.5 

3.1 

2.9 

^ < *t 

it 

3 

0.119 

0.144 

22.6 

11.0 

12.7 

4.9 

4.8 

5.2 

" 

5 

0.092 

0.216 

26.6 

20.7 

15.6 

6.4 

5.4 

7.2 

If 

6 

0.070 

0.128 

25.2 

12.4 

5 . 5 

3.9 

5.1 

2.7 

'* 

7 

0.059 

0.112 

15.2 

6.5 

4.0 

2.0 

3.0 

2.0 

*' 

8 

0.108 

0.185 

31.9 

10.6 

5 . 6 

3.8 

1.9 

2.0 

76174 

1 

0.055 

0.130 

16.6 

15.3 

14.4 

3.7 

3.7 

3.0 

'* 

2 

0.234 

0.257 

32.5 

11.9 

8.4 

5.4 

3.8 

3.3 

76175 

2 

0.103 

0.187 

46.2 

35.3 

29.6 

3.9 

8.8 

4.9 

'* 

3 

0.038 

0.102 

23.4 

17.8 

19.0 

4.3 

6.1 

4 . 5 

4 

0.151 

0.131 

54.2 

23.4 

27.3 

11.1 

2.6 

3.5 

6 

0.128 

0.152 

28.4 

10.9 

10.7 

4.9 

6.3 

3 . 6 

** 

7 

0.191 

0.223 

57.1 

31.9 

32.3 

2.6 

1.8 

2.  3 

II 

3 

0.124 

0.187 

45.5 

13.8 

3.7 

4.1 

1.4 

2.0 

*' 

9 

0.151 

0.245 

49. 8 

20.0 

4.3 

3.6 

3.6 

2.0 

76197 

1 

0.046 

0.074 

8.8 

3.6 

2.8 

2.7 

2.0 

1.3 

75338 

1 

0.125 

0.304 

32.9 

21.5 

9.7 

3.7 

3.4 

3.4 

II 

2 

0.070 

0.105 

15.1 

2.8 

4.3 

2.8 

1.4 

1.  8 

" 

3 

0.115 

0.277 

33.1 

26.1 

20.6 

4.6 

5.9 

2.7 

4 

0.094 

0.171 

21.9 

10.4 

10.0 

14.3 

14.0 

10.3 

76162 

1 

0.160 

0.303 

39.5 

17.2 

15.5 

9.5 

13.8 

9.0 

'* 

2 

0.165 

0.  194 

45.2 

12.2 

8.0 

17.3 

3.9 

3.1 

3 

0. 172 

0.226 

43.3 

7.9 

5.6 

5.  3 

2.7 

2.1 

" 

4 

0.157 

0.287 

34.7 

8.9 

10.4 

14.5 

8.1 

7.7 

II 

5 

0.218 

0.259 

68.0 

7.3 

9.9 

11.6 

2.9 

2.9 

6 

0.159 

0.317 

58.4 

20.1 

5.9 

15.0 

3.9 

2.3 

It 

7 

0.150 

0.258 

41.3 

9.0 

6.6 

14.2 

6.2 

5.2 

106 


li 

that  the  Sight  Eval  data  does  indeed  present  a substantially  true 
picture  of  the  system  performance,  the  data  conveys  a strong 
impression  of  erratic  system  performance. 

It  seems  reasonable  to  conclude  that  whatever  the  causes 
of  the  radar  lag  problem,  the  lag  angle  is  not  a function  of  the 
engagement  geometry.  This  cor.clusion  can  be  criticized  on  the 
basis  that  a higher  order  fit  should  have  been  tried  or  some 
conceptually  different  type  of  fit  involving  different  engagement 
variables  might  have  been  successful,  etc.  Perhaps.  But  the 
fact  remains  that  the  various  turn  rates  used  all  have  roughly 
the  same  direction  as  the  vector  lag  angle  and  most  of  the 
engagement  variables  are  involved  in  computing  these  turn  rates. 

There  should  have  been  some  perceivable  trend  in  the  fit  constants 
somewhere.  After  all,  the  lag  angle  problem  is  caused  largely 
by  an  incorrect  axis  turn  rate  with  perhaps  an  initial  misorientation 
problem  added  for  good  measure. 

There  does  exist  the  possibility  of  taking  one  of  the  fits 
and  working  up  those  values  of  its  constants  which  yield  the  best 
overall  statistical  decrease  in  radar  lag  angle  for  a large 
number  of  engagements.  But  there  are  two  obvious  hazards  in  the 
procedure.  First,  if  the  true  source  of  the  problem  remains 
undiscovered,  what  assurance  is  there  that  there  will  not  be 
radical  changes  in  system  performance  in  the  future?  Second, 
since  the  engagements  are  so  similar  and  so  brief,  would  the 
correction  be  valid  in  engagements  having  a different  geometry 
or  a greater  duration? 

One  further  comment  is  necessary.  It  may  seem  that  the 
more  direct  approach  of  doing  a system  analysis,  on  the  MA-1 
system  and  its  subsequent  modifications,  was  ignored.  After  all, 
the  original  applications  of  the  MA-1  system  radar  did  not  require 
it  to  produce  the  highly  accurate  target  angular  direction  data 
now  being  asked  of  it.  An  analysis  of  the  MA-1  system  design 
and  performance  was  initiated.  This  analysis  required  detailed 
system  design  theory  and  performance  data  as  well  as  systems 
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modification  and  modified  system  performance  data  of  the  MA-1 


used  in  Sight  Eval  testing.  The  extensive  period  since  the 
original  design  and  test  of  the  MA-1  materially  reduced  the 
availability  of  pertinent  design  and  performance  data,  leaving 
the  geometrical  approach  as  the  only  analysis  alternative. 

1 

1 


3. 3. 4. 3 Time  Data  Phase  Shift 


r 


An  examination  of  the  plots  of  radar  antenna 
position  in  azimuth  and  elevation  versus  film  target  position 
shows  an  apparent  data  phase  shift  which  varies  up  to  .5  seconds 
and  averages  approximately  .25  seconds.  There  are  a number  of 
reasons  why  such  a phase  shift  could  exist.  Delay  in  instru- 
mentation circuitry  is  a function  of  the  number  of  processing 
and  amplification  stages  through  which  an  electronic  signal 
must  pass.  However,  this  category  of  delay  results  in  less 
than  a microsecond  per  stage  and  will  seldom  reach  the  .25 
second  delay  observed.  A more  likely  possibility  is  a discrepancy 
in  the  synchronization  of  magnetic  tape  and  film  data  during  the 
merger  into  a single  data  tape.  Of  course,  the  probability  is 
high  that  the  phase  shift  does  indeed  exist,  and,  if  so,  the 
direct  modification  to  the  correction  algorithm  is  obvious. 

Graphs  of  passes  that  illustrate  the  data  shift  are  shown  in 
Figures  33  thru  37. 

The  most  obvious  example  of  a data  shift  is 
shown  in  Figure  33.  On  this  pass  the  radar  was  leading  the  film 
data  in  azimuth.  This  immediately  suggests  a data  merging  or 
calibration  problem  as  the  radar  can  never  constantly  lead  the 
target.  The  occurrence  of  radar  lead  was  not  often  obvious, 
but  Figure  34  strongly  suggests  that  the  radar  was  leading 
the  elevation  for  approximately  1.1/2  of  the  3 seconds  of 
available  data.  The  radar  data  follows  the  film  data  profile 
very  closely  and  the  differences  between  film  and  radar  are 
small,  and  rarely  exceed  10  MRAD.  If  not  closely  inspected, 
this  pass  would  have  been  accepted  as  an  ideal  pass. 

The  third  example  of  a possible  data  shift  in 
azimuth  is  seen  in  Figure  35.  The  elevation  radar  data  closely 
follows  the  film  data  wxth  the  radar  occasionally  overshooting 
the  target,  but  immediately  being  corrected.  This  is  how  one 
would  expect  the  radar  to  behave,  and  the  component  differences 
rarely  exceed  15  MRAD.  Inspection  of  the  azimuth  data  reveals  a 
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Figure  35  Lag  in  Radar  Azimuth 


very  poor  tracking  job  by  the  radar.  The  same  basic  profile  is 
observed  by  radar  and  film  but  the  radar  is  lagging  the  film 
by  approximately  25-30  MRAD.  Also  the  radar  profile,  if  smoothed, 
seems  to  be  about  1/4  second  out  of  synchronization  with  the 
film. 

The  above  examples  are  extremes  noted  in 
reviewing  the  data.  The  largest  percentage  of  the  passes  showed 
the  radar  lagging  the  film  data.  This  lag  time  was  felt  to  be 
excessive  in  most  cases.  The  radar  profiles  in  Figures  36  and  37 
suggest  the  radar  is  roughly  1/4  second  behind  the  film  data 
and  are  typical  of  many  passes.  Being  unable  to  isolate  the 
true  radar  lag  from  a data  shift,  it  was  decided  to  try  a 
variety  of  time  shifts.  It  was  felt  that  the  radar  should 
respond  with  a correction  in  less  than  a tenth  of  a second. 
Accordingly,  time  corrections  of  4 to  8 frames  were  selected 
to  see  if  the  total  lag  angle  could  be  consistently  corrected. 

Program  SDATA  described  in  detail  in  Section 
3. 3. 4. 2 was  modified  to  perform  all  lag  angle  calculations  in 
the  program  inertial  coordinate  system,  with  the  film  target 
direction  being  transformed  to  inertial  coordinates  by  that 
transformation  consistent  with  the  lag  time  assumed.  Two  tapes 
were  run  with  this  correction.  The  results  of  the  data  shift 
are  shown  in  Tables  14  and  15. 

In  many  passes  the  lag  angle  is  reduced  below 
the  15  MRAD  requirement.  In  Table  14,  5 of  the  29  passes  had 
lag  angles  less  than  15  MRAD.  After  the  frame  shift,  24  of  29 
had  lag  angles  less  than  15  MRAD.  In  Table  15,  4 of  the  26  passes 
had  lag  angles  less  than  15  MRAD.  After  the  frame  shift,  14  of 
24  had  lag  angles  less  than  15  MRAD.  Of  the  54  passes  tried 
for  this  data  shift,  only  2 that  would  have  been  accepted  were 
made  unacceptable  by  the  shift.  However,  there  appears  to  be 
a definite  contribution  to  lag  angle  from  this  data  shift. 

Further,  an  error  in  synchronization  is  the  only  possible  explana- 
tion for  those  few  passes  which  show  the  radar  leading  the  film. 
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TABLE  14  Comparison  of  Average  Lag  Before  and  After  Frame  Shift 


BEST 


AVERAGE 

FRAME 

ADJUSTED 

FRAMES 

LAG 

SHIFT 

LAG 

ERROR 

MISSION 

PASS 

MRAD  (No 

of  Frames)  MRAD 

< 15  MRAD 

COMMENTS 

76111 

1 

23.2 

4 

6.12 

4-6 

2 

21.8 

6 

7.73 

4-8 

3 

40.6 

5 

4.50 

4-7 

4 

36.8 

5 

4.34 

4-6 

5 

45.8 

6 

2.07 

5-7 

76118 

1 

28.8 

6 

3.93 

4-8 

2 

30.8 

5 

6.65 

4-6 

3 

42.0 

6 

5.44 

5-8 

4 

15.4 

4 

14.89 

4 

5 

33.6 

5 

7.23 

4-6 

6 

74.9 

6 

5.97 

5-6 

7 

66.5 

6 

15.03 

— 

8 

135.5 

8 

130 .48 

— 

No  Radar  Lock 

76119 

1 

383.5 

4 

388.14 



No  Radar  Lock 

2 

11.8 

8 

4.02 

4-8 

3 

16.5 

4 

6.18 

4-7 

4 

25.8 

6 

7.96 

4-8 

5 

20.6 

6 

16.81 

— 

6 

36.4 

8 

12.85 

7-8 

7 

45.9 

6 

5.59 

5-8 

76121 

1 

5.7 

4 

7.28 

4-8 

2 

11.8 

4 

10.11 

4-8 

3 

10.7 

6 

6.79 

4-8 

4 

14.6 

4 

9.99 

4-6 

5 

50.3 

6 

6.71 

5-8 

6 

32.3 

4 

11.29 

4 

7 

24.1 

4 

3.92 

4-6 

8 

17.6 

4 

19.33 

— 

9 

55.1 

6 

4.57 

00 

1 

m 

TABLE  15 


Comparison  of  Average  Lag  Before  and  After  Frame  Shift 


BEST 


AVERAGE 

FRAME 

ADJUSTED 

FRAMES 

LAG 

SHIFT 

LAG 

ERROR 

MISSION 

PASS 

MRAD  (No 

of  Frames)  MRAD 

< 15  MRAD 

COMMENTS 

76077 

1 

10.7 

4 

6.19 

4-8 

2 

14.9 

4 

6.34 

4-8 

3 

21.6 

5 

10.33 

4-6 

4 

8.3 

4 

18.75 

— 

5 

30.2 

5 

7.88 

4-6 

6 

15.9 

4 

25.61 

— 

76078 

1 

8.3 

4 

6.13 

4-8 

2 

23.4 

5 

6.25 

4-8 

3 

15.9 

4 

14.73 

4-6 

4 

17.5 

8 

8.19 

4-8 

5 

18.8 

4 

16.02 

— 

6 

47.5 

7 

8.99 

5-8 

7 

35.9 

5 

11.92 

4-6 

8 

45.0 

6 

8.07 

4-7 

76083 

1 

37.5 

5 

8.81 

4-7 

2 

30.2 

4 

16.64 

— 

3 

33.7 

4 

6.45 

4-5 

4 

37.6 

5 

13.88 

4-5 

76085 

1 

26.5 

7 

20.78 

— 

76092 

1 

35.0 

4 

34.22 

• w 

2 

13.01 

4 

24.07 

— 

Radar  Lead 

(Elevation) 

76093 

1 

25.1 

8 

15.52 

— 

2 

63.5 

4 

78.92 

— 

No  Radar  LocJc 

3 

185.3 

8 

181.59 

— 

No  Radar  Loc)c 

4 

226.6 

8 

225.84 

— 

No  Radar  Lock 
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Again  reaching  a general  conclusion  is  made 
difficult  by  the  lack  of  detailed  documentation  on  the  experi- 
mental design  and  on  the  data  reduction  procedures.  This  pro- 
cedure did  not  give  a great  improvement  in  all  passes,  but  the 
results  do  indicate  that  the  radar  may  be  considerably  more 
accurate  in  determining  target  angular  direction,  than  was 
previously  thought. 

3. 3. 4. 4 Suggestions  for  Future  Study 

It  is  evident  that  a time  lag  occurred  between 
the  times  the  film  and  radar  data  were  recorded.  This  lag  may 
be  an  actual  interval  introduction  by  real  radar  error  or  it 
may  be  an  artificial  interval  caused  by  instrumentation  or  data 
merger  into  a single  tape.  This  time  lag  which  results  in 
radar  error  is  of  sufficient  importance  that  a follow-on  study 
should  determine  the  true  source  of  the  delay. 

In  connection  with  the  current  analysis,  the 
availability  of  a detailed  instrumentation  and  test  plan, 
maintenance,  calibration,  and  modification  logs  for  the  aircraft 
and  fire  control  system,  and  details  of  data  processing  would 
have  been  of  material  assistance.  In  future  testing,  the  close 
coordination  between  analysis  and  test  planning  should  be 
emphasized,  with  consideration  given  to  a single  analysis/Test 
Manager.  A sample  analysis  exercise  concurrently  with  initial 
flight  tests  will  eliminate  many  later  analysis  snags. 
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APPENDIX  B 


PROGRAMS  FOR  KALMAN  FILTERS 


PART  A - CONSTANT  ACCELERATION  FILTER 


SECTION  PAGE  NO. 

DESCRIPTION 

METHOD 

INPUT 

OUTPUT 

PROGRAM  LISTING 

PART  B - CONSTANT  VELOCITY  FILTER 

DESCRIPTION 

METHOD 

INPUT 

OUTPUT 

PROGRAM  LISTING* 

*Listing  is  of  the  Straight  Line  Trajectory  Model 
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DESCRIPTION 


Computer  program  'KALMAN'  is  a FORTRAN-coded  range 
type  Kalman  filter.  It  contains  the  necessary  equations  to 
use  one  of  three  models;  constant  acceleration,  sinusoidal 
acceleration  (2G  amplitude)  and  sinusoidal  acceleration 
(6G  amplitude) . In  addition,  the  Kalman  gain  may  be  computed 
in  one  of  three  ways;  time  varying  gain  (gain  computed  each 
pass  as  function  of  the  filter  parameter ' P ') , synthetic  time 
varying  gain  (gain  computed  each  pass  as  function  of  time 
and  input  coefficients) , and  constant  gain  (gain  is  input) . 

The  program  computes  statistics  which  can  be  used  to  com- 
pare range,  velocity  and  acceleration  with  their  estimated 
values . 

The  purpose  of  the  model  is  two-fold: 

1)  To  provide  a means  of  evaluating  several  variations 

of  the  filter. 

2)  To  enable  the  user  to  determine  optimum  values 

of  some  significant  filter  parameters 

The  program  is  not  and  never  was  intended  to  be  the  end 
product  of  the  filter  study. 


154 


r 


METHOD 


I 


Calculate  'True'  Values  of  Range,  Velocity,  and  Acceleration 
(R,  R,  R)  According  to  One  of  Three  Methods: 

k is  an  integer  which  takes  on  all  values  from 
1 through  NPTS  (NPTS  = duration  of  run/sampling  rate 

T = sampling  rate,  sec. 

a - Constant  Acceleration  Model 

R(k)=R^+V^(kT)+  I A(kT)^ 

R(k)=V  +A(kT) 
o 

R(k)=A 

b - Sinusoidal  Acceleration,  2G  Amplitude 

R(k)-R^+24.9  55[ (kT) -1.27 324SIN(. 0 31416  k)  ] 
R(k)=24.955[1.0-COS(. 031416  k) ] 

R(k)=19.60SIN(. 031416  k) 
c - Sinusoidal  Acceleration,  6G  Amplitude 

R(k)=R^+2.99466 [k- 31 . 831SIN (o. 31416  k) ] 
R(k)=74.8665[1.0-COS(. 031416  k) ] 

R(k)=58.8SIN(. 031416  k) 


Generation  of  Gaussian  Noise: 

Generate  12  random  numbers  each  having  a 
uniform  distribution  between  0 and  1 

(represented  by  RN . ) 

1 12 

Then  the  noise,  n(k)=  ®v  [ ^RN  . ) -6 . 0 ] is  a 

i=l 

Quasi-Gaussian  random  variable  with  zero  mean 
and  standard  deviation=cJv 

The  'Measured'  Value  of  Range  is  Z (k) =R (k) +n (k) 


(4) 


(5) 


X.  , X„,  X^=  Estimates  of  R,  R,  R 


2^  3 

To  initialize: 


X^=Z(3) 


X2  = 12.5[3.0Z(3)-4.0Z{2)  + (1)  ] 


X2=0. 0 


Filter  Equations 

a)  P(k|k-1)=  ((J))  P(k-l|k-l)  ((1)T)+S(k) 

^ 1 T T^2.0 

|)=  Constant  Matrix  = 


1 

0 


T 

1 


(()T=  Transpose  of  ({) 

^w^T^o  2 

S(k)=  Maneuver  Noise  Term=  / v 


P is  initialized  as 


a 2 

V 


(1/T^) 


b)  Gain  (Time  Varying  Gain) 


K(k)  = 


1.0 


Pj^^  j^(k|k-l)+''v 


P^^^(klk-l) 

P2\{k|k-1) 

P^'j^Cklk-l) 


c)  Estimates: 


3T- 


6T 


9T 

18T 


6T 

18T 

36,0 


1 (3/2T) 

(3/2T)  (13/2T^) 


(6/T^) 


(1/T^) 

(6/T^) 

(6/t'^)+°w^ 


x(k|k)  = ((j))  x{k  |k-l)+K(k)  [Z(k)-(H)  ( (j))  X (k  |k-l)  ] 
where  (H)  = Row  Matrix  =(1,0,0) 


d)  P(k|k)=[T  - K(k) (H) ] P(k|k-1) 


where  T = Identity  Matrix  = 


10  0 
0 10 
0 0 1 
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Note  - For  synthetic  time  varying  gain  and  constant 
gain  equations  (a)  and  (d)  are  By-Passed  for 
the  constant  gain  case,  Ki=input  values  for 
synthetic  time  varying  gain,  ••• 

C7 

+ C6;  K3  = 


= 


Cl 


+ C3;  K„  = 


C4 


k+C2  ^ k2+C5  k3+C8 

where  Cl,  C2, ,C9  are  input 


+ C9 


(6)  Statistics  - The  Following  RMS  Deviations  are  Calculated: 


/ 1.0 

N2 

E 

[R(k)=X3(k| 

|k)]2 

^ (N2-N1+1) 

k=Nl 

1 1.0 

N2 

E 

[R(k)-X2 (k] 

|k) 

V 

1 (N2-N1+1) 

k=Nl 

(AX2)RMS  = 


1.0 


N2 

E 


(N2-N1+1)  k=Nl 


[R(k)-X3(k|k) 


Nl,  N2  are  allowed  to  take  on  3 sets  of 
values,  one  of  which  is  always  Nl=4; 
N2=NPTS.  The  other  two  sets  are  input. 


!,  \ 


i ‘ 


i i 


■ i 


i 


PREPARATION  OF  INPUT  CARDS 

Note  •••  F.PT.  = Floating  Point,  I=Integer 
CARD  TYPE  #1  (Mandatory) 


Var.  MNEMONIC  COLS.  TYPE 
IMEAS  5 I 


DEFINITION 
Model  Indicator: 

0 = Constant  Acceleration, 

R,  R,  R computed 

1 = Constant  Acceleration, 

R Measurements  Input 

2 = Sinusoidal  Acceleration, 

2g  Amplitude 

3 = Sinusoidal  Acceleration, 

6g  Amplitude 


IGAIN 

10 

I 

Gain  Indicator; 

0 = Time  Varying  Gain 

1 = Constant  Gain 

2 = Synthetic  Time  Varying 

Gain 

INTIA 

11-15 

I 

Statistics  Interval  #1,  First 

Point 

INTIB 

16-20 

I 

Statistics  Interval  #1,  Last 

Point 

INT2A 

21-25 

I 

Statistics  Interval  #2,  First 

Point 

INT2B 

26-30 

I 

Statistics  Interval  #2,  Last 

Point 

■^R 

TRATE 

31-40 

F.PT. 

Sampling  Rate,  Sec. 

SECS 

41-50 

F.PT. 

Total  Time  of  Run,  Sec. 

(SECS^^j^^g  _ 250  Maximum) 

CARD 

TYPE  #2 

(Mandatory) 

0 2 
V 

SVSQ 

1-10 

F.PT. 

Variance,  Measurement  Noise 

a 2 
w 

SWSQ 

11-20 

F.PT. 

Variance,  Manuever  Noise 

r 

RO 

21-30 

F.PT. 

Initial  Range,  Meters 

o 

(Not  required  if  IMEAS=1) 

V 

VO 

31-40 

F.PT. 

Initial  Velocity,  M/SEC. 

(Not  required  if  IMEAS=2,3) 
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1 

a 

A 

41-50 

F.PT. 

2 I 

Acceleration,  M/SEC  | 

(Not  required  if  IMEAS=2,3)  j 

' ■ 

51-60 

F.PT. 

Kalman  Gain  (Required  for  ; 

IGAIN=1  Only)  1 

^2 

^2 

61-70 

F.PT. 

Kalman  Gain  (Required  for 

IGAIN=1  Only)  ' I 

^3 

^3 

71-80 

F.PT. 

Kalman  Gain  (Required  for  ! 

IGAIN=1  Only)  ! 

j 

CARD 

TYPE  #3  (This 

Card  Required 

for  Synthetic  Time  Varying  Gain  I 

Only 

IGAIN=2) 

i 

i 

Var. 

MNEMONIC 

COLS. 

TYPE 

DEFINITION  1 

i 

Cl 

1-10 

F.PT. 

Coefficients See  Equations  | 

On  Page  3 1 

C2 

11-20 

F.PT. 

i 

^3 

C3 

21-30 

F.PT. 

^4 

C4 

31-40 

F.PT. 

C5 

41-50 

F.PT. 

^6 

C6 

C7 

51-60 

61-70 

F.PT. 

F.PT. 

i 

1 

^8 

C8 

71-80 

F.PT. 

1 

I 

S 

C9 

1-10, 
2^^  CD. 

F.PT. 

} 

J 

5 

I 

j 

CARD 

TYPE  #4  (This 

Card  Required 

I 

for  IMEAS=1  Only)  1 

i 

Enter  successive  values  of  range  measurements  8 per  card, 

8F10.1  Format.  Use  number  of  cards  required,  but  number 
of  values  must  not  exceed  250. 
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The  program  output  consists  of  three  (3)  parts: 

1)  Headings  (type  of  run,  values  of  input  parameters, 

..AAA 

etc.).  Plus  the  values  of  R,  R,  R,  X^,  X^,  K^, 

at  each  point. 

2)  The  measurement  and  noise  at  each  point. 

3)  Statistics  computed  at  three  intervals,  two 
of  which  are  input  and  the  last  of  which  is 
for  the  range  (4  — ♦■NPTS)  . 


DESCRIPTION 


r 


Computer  program  'CNVEL'  is  a FORTRAN  coded  range  type 
constant  velocity  Kalman  Filter.  The  program  contains  the 
necessary  equations  to  implement  one  of  three  models  (trajectories) : j 

constant  accelerations,  straight  line,  sinusoidal  acceleration  ■ 

(2g  amplitude) , and  a sinusoidal  acceleration  (6g  amplitude) . 

The  program  computes  statistics  which  can  be  used  to  compare 
range,  velocity  and  acceleration  with  their  respective  estimated 
values . 

The  purpose  of  the  program  is  twofold  - 

1)  to  provide  a means  of  evaluating  several 
variables  of  the  filter, 

2)  to  enable  the  user  to  determine  optimum 
values  of  some  significant  filter  parameters. 
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L. 


METHOD 


(1)  Calculate  'true'  values  of  range  velocity  and  acceleration 
(R,  R,  R)  according  to  one  of  three  methods: 


b . 


c . 


Constant  acceleration,  : 

R^  = Xj^(k) 

= X^(k-l)  + ■ 

R^  = (k) 

= X2 (k-l)  + ' 

R.  = X3(k) 

= constant 

Sinusoidal 

acceleration 

R^  = X^(k) 

= 0.9982  tk 

R.  = X^<k> 

= 24.955  [1. 

s’i  - 

= 19.60  sin 

Sinusoidal 

acceleration 

R^  = X^(k) 

= 2.9947  [k 

Ri  = Xj  (k) 

= 74.866  [1. 

Rj  = Xj(k) 

= 58.8  sin  ( 

T 

2 


where  T = sampling  rate,  sec.,  and 

k = is  an  integer  which  takes  on 
all  values  1 through  NPTS 

(NETS  = T_  - T )/At 
r O 


(2)  Generation  of  Gaussian  noise: 

Generate  12  random  numbers  each  having  a uniform  distri- 
bution between  0 and  1 - represented  by  RN . 

12 


Then  - noise  = n(k)  =a  ZRN.  -6.0i 

V L 1 J 


.s  a quasi-gaussian 
random  variable  with  zero  mean  and  standard  deviation  = a . 


(3)  The  range  measurement  is  Z (k)  = R(k)  + n(k) 


167 


(4)  X^,  X^,  = estimates  of  R,  R,  R 

To  initialize: 


= Z(2) 

X2  = 25.0  [Z  (2)  - X(l) ] 

X^  = 0.0 


(5)  Filter  Equations 

a.  Initialization 

P(2/2)  = 625  a. 


b.  Iteration 

P(k  + l|k)  = (f)P(k|k)  4)'+  Q' 
where : 


2 

. 0016 

.04 

.04 

2.0 

4,'  = 


Q'  = a 


1 

0 

1 

,04 


w 


. 04 
1 
0 
1 


t'^/4 

T^/2 


T^/2 


K(k  + 1)  = 


P^^{k  + ilk)  + 


Pj^^(k  + l|k) 

P23^(k  + ilk) 


X(k  + ilk  + 1)  = 4>  X(klk)  + K(k  + 1)  Z (k  + 1)  - {1.04)X(klk) 


P(k  + l|k  + 1)  = 


1-K^(k  +1)  0 


-K2(k  +1)  1 


P(k  + ilk) 
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The  following  RMS  deviations  are  computed: 


i 


where  N1  7^  N2  = initial  and  final  indices  of  the 

points  to  be  utilized  in  the  computations 
of  statistics 


— — 

i 

INPUT  PREPARATION 


CARD  TYPE  #1 


VARIABLE 

COLS. 

TYPE 

DEFINITION 

R 

o 

1-10 

F.Pt. 

Initial  range,  meters 

R 

o 

11-20 

F.Pt. 

Initial  velocity,  meters/sec 

R* 

o 

21-30 

F.Pt. 

2 

Initial  acceleration,  m/sec 

T 

o 

31-40 

F.Pt. 

Initial  time,  sec 

41-50 

F.Pt. 

Final  time,  sec 

51-60 

F.Pt. 

Time  increment,  sec 

CARD  TYPE 

#2 

1-10 

F.Pt. 

Variance,  maneuver  noise 

11-20 

F.Pt. 

Variance,  measurement  noise 

NEND 

21-25 

Int. 

Value  of  index  of  highest 
point  to  be  printed 

CARD  TYPE 

#3 

NOST 

1-5 

Int. 

Number  of  sets  of  statistics 
to  be  computed  (3  in  most  cases 

CARD  TYPE 

#4 

(1  card 

per  statistics  set) 

NST 

1-5 

Int. 

Initial  and  final  'K'  indices 
of  points  to  be  used  in 

NFN 

6-10 

Int. 

computation  of  statistics  set 

■j 

\ 
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OUTPUT 


I 


2 2 

(1)  Input  values  of  R , R , R , o , o 

O O O V w 

(2)  Iterative  output  at  each  'k' 

k,  R(k)  , R(k)  , R(k)  , Xj^(k),  X2(k),  ^^(k), 
xkA,  xkB  (Kalman  gaon  constants) , 
z(k)=  measurement  of  range 
xN(k)=  noise 

(3)  Statistics:  RMS  deviations  for  position  and  velocity 


I 


fi 


:i 


u- 

5 ..  . UJ-< 


CT  r'j 

» • ■7' ' 


aLn  ^itj  O'^ 
— 'C  ti»-<*— 
iQCiiS— Cl 


•{  >i4  ^ O ^ 

-'xu*Mh-ir 


w3  ,-5  C ^ 

'C  <e  ;t  j oj 

o U-  X X -13  it 

cti  k-i  M X 'j  :• 


•cocv^j  or 
c 'tn  w o Ci 

2S  I)  :«t3  r^  r 

♦—  t^O  'J  » 

— 0 — ' ‘l  • • 

■f_(  W r. 

-t;CJ  < t-t 

C iH 

XM  «a  II 

II  a^u 


05  r:  *■ 
X 3 

30  *04 

3 6-  e?-'- 

-TSCU  .:bh 

-X3>ace- 

C II  *5  II  l|  C! 

s O 5-OM  ^ 

I 

f5 » O — a*  - 

5XUo:x^ 


6-6-23 
CUCl3» 
a XO' 
: ♦ ♦ • 

3 OCC 
c ■-  • •in 
- CLjininr>4 
-tXfNrs^' 
S II  II  11  tl 

3 .e  aua 

; Q«Q4C^a« 


Qr-i  ftiOO 

Ca*-  O^OO 

•J  C-croo 

OOC. 

<•  o • • 

L3COOJ 

saxu*!  4 '«  •«>  < 

tu  % <•  ..«^oa» 

aB3isi**>  ^a^a^Uiia 

--J-4  ♦ ♦ ♦ ^ 

f3  •/— ' % 

ooooo 

w^ucu  « « . . 


'tf  < a ra  I 

XXHXX 


3 OOH  *1^—53  •C'CCC 

I M -- «-r-  * *iC  ♦ 

■— oo  atfaj<H-«-<auQ 
• m I 6-  d»5 11  II  ^ 1C  6-  a«  a*  X 
-~-Tc I-* « X < a t-j-sx  u n II  II 


♦ M 

oA«ca 

6-^  •^C'O^  **-  ■- 

WW  ox*  *00  IT  V3 

xa  4 ♦ «..^oio  o o 

4 4 < 4 ^ n R 

•<<  2i  « 

OC*  ■^'O'  oaT^ 

— — ^X|  |»*CC^Mll 

NV.  xx^oaxa 
«:a  M It  • t^xaiTH 
Oo**-  M 

H II  4*-J-4'-~-'|  I VQlCX- 
•OK— II  II  II  i|  st<a 
KK  II  < aOUQOtOOO 
XX-JXXD«CUC-a4!J  X6uO 


_ — - ' 

i 

t 

1 

The  attached  dociamentation  pertains  to  the  program 
'ERRORV  (velocity) . The  decks  for  two  additional 
i programs;  'ERR0R'  (acceleration)  and  'ERR0RR' 

(range)  are  also  attached.  The  three  programs  are 
; very  similar.  The  following  is  a list  of  the  distri- 

i butions  utilized  in  each  of  the  three  programs  - 


PROGRAM 

VARIABLES  TYPE  DISTRIBUTION 

'ERR0RV' 

^t'  ^a 

uniform 

range:  200,  600  KTS 

V.  &v  errors 
t a 

normal 

a = 4.0  KTS,  x=0 

“t'  “a 

uniform 

range:  -5°,  +20° 

a &a  f v"‘ors 
u a 

normal 

Q 

It 

O 

M 

0 

o 

roll  angle ((j)) 

uniform 

range:  -90°,  -60°  & +60°, 

pitch  angle (0) 

uniform 

range:  -30°,  +45°  (75%) 

-90°,  -30°  (12.5%) 
+45°,  +90°  (12.5%) 

yaw  angle  (i))) 

uni from 

range:  -180°,  +180° 

sine,  cosine 
errors 

normal 

a = 0.005,  x=0 

' ERROR ' 

at  , aa 
x'  X 

uniform 

a = 32.16,  x=0 

at  , aa 

y y 

uniform 

Q 

II 

00 

• 

o 

o 

a-comp  errors 

normal 

pdf=-8. 0+4e  ^^^^^where 
x=normal  dist,  o=1.0&x=0 

(same  as 

above) 

sine,  cosine 
errors 

(same  as 

above) 

'ERR0RR' 

e ,n 

uniform 

range:  -15°,  +15° 

e,ri  errors 

normal 

0 = .0691°,  x=0 

range 

uniform 

range:  0,  2500  ft. 

range  error 

normal 

0 = 4.0  ft . , x=0 

1 


+90® 
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COMPUTER  PROGRAM  ' ERRORV ' CONTAINS  THE 
MONTE  CARLO  CALCULATION  OF  THE  ERROR 
DISTRIBUTION  OF  THE  VELOCITY  OF  THE 
TARGET  RELATIVE  TO  THE  ATTACKER  IN  THE 
ATTACKERS  BODY  FRAME  OF  REFERENCE. 


Page  2 

Pages  3,  4,  5 


Pages  6,7,8 


Description  of  lnput/0utput 

Explanation  of  program  cross-referenced  to 
line  numbers  on  compilation  listing 

Purpose  and  method  used  in  each  subroutine 


The  Program  Requires  56000  Words  (Octal) 
of  Storage  on  the  CDC  6600. 


COMMENTS 


LINE  NOS. 

7-11 

12 

15 

16-18 


20-21 

23 

25 

27-29 

31-33 

35-37 

39-41 

43-48 


50-52 

53-59 


Define  Constants 

READ:  Seed  for  random  number  generator 

and  number  of  passes 

Test  for  end  of  file  on  Input  Tape  5 

If  last  case  has  been  processed,  print 
'NORMAL  TERMINATION'  and  stop -if  another 
case  is  to  be  processed,  proceed  to 
statement  No.  30 

Print  Input  Values 

Initialize  Random  Number  Generator 

Start  Main  Execution  Loop 

Select  airspeeds  v,  & v from  a uniform 
distribution  having  the^range  200-600  KTS 

Select  airspeed  erro.vs  from  a normal 
distribution  having  a a=4.0  KTS  and  a 
MEAN=0 

Select  angles  of  attack  a & from  a 
uniform  distribution  having  a range 
-5°  to  +20°  (convert  to  radians) 

Select  angle  of  attack  errors  from  a 
normal  distribution  having  a a=0.1° 
and  MEAN=0.  (Convert  to  radians) 

Compute  airspeed  components,  target 
and  attacker 


Select  yaw  angles  from  a uniform  distri- 
bution with  the  range  -180°  to  180° 
(convert  to  radians) 

Select  roll  angles  from  a uniform  distri- 
bution with  the  range  =90°  to  -60°  and 
+60°  to  +90°  (convert  to  radians) 
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LINE  NOS. 


COMMENTS 


60-71 


73-80 

82-84 

85 

91 

93-95 

96 

97-99 

100 

101-110 


112 

114-115 

116 

117-119 


120 

122-125 

126 


Select  pitch  angles  from  uniform  distri- 
bution with  the  range:  -90°  to  +90° 
where 


-90°  to  -30°  Occurs 
-30°  to  +45° 

+45°  to  +90° 


12  1/2  % of  the  time 
75%  **  **  ** 

12  1/2% 


(convert  to  radians) 

Select  sine  & cosine  errors  from  a normal 
distribution  with  a=  .005,  MEAN=0. 

Set-up  'Y'  matrix  terms 

Compute  elements  of  'Y'  matrix 

Solve  matrix  equation,  no  error  terms 

Set-up  'Y'  matrix  terms  (errors  included) 

Compute  elements  of  'Y'  matrix 

Set-up  'X'  matrix  terms  (error  included) 

Compute  elements  of  'X'  matrix 

Add  error  terms : 


< 

rt 

= v^  + e 

7 

V 

a 

= V + e 
a 

“t 

= + e 

1 

1 

a 

a 

= a + e 
a 

^tx’ 

= cosa^ 

^t 

f 

V 

ax 

= cosa 

a 

^ty' 

1 

= 0 

1 

1 

r 

''a/ 

1 

= 0 

^tz 

= sina^ 

/ 

V 

az 

= sina 

cL 

V 


Solve  matrix  equation  (with  error  terms) 


V =[X'] 
V 


[Y'l 


[Vt'l  - [v^'l 


Difference  velocity  components  (v^'  - v^) 

Compute  total  velocity  defference  from 
sum  of  squares 


Store  (v 


vx 


- V ) 
vx 


(Vvy’  - Vvy)  , (Vvz'-  Vvz) 


and  output  array 

End  Main  Execution  Loop 
Print  Output 

Return  to  statement  No.  1 for  next  case 
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SUBROUTINE  'RANDOM' 


Purpose:  To  compute  'NRN'  random  numbers  from  either 

a uniform  or  normal  distribution. 

Subroutine  'RANF'  is  utilized  to  generate 
a uniformly  distributed  random  number  (0-»l) 


IDIST=0  Uniform  Distribution 
Cl  = Additive  Term  (To  Scale  Range) 

C2  = Multiplier  Term  (To  Scale  Range) 


IDIST=1  Normal  Distribution 
Cl  = 2.  - STD  Deviation 

C2  = X = Mean 


Method: 

a)  Uniform  Distribution 


Start  Loop  'A' 
i=l,  NRN 

I 

RN=RANF (DUM) 

i 

RNS.=  (RN+Cl)  C2 

^ i 

End  Loop  'A' 


b)  Normal  Distribution 


Start  Loop  'B' 
i=l,  NRN 

i 

RNT0T=O 
Start  Loop  'C 

j=l,  12 

RN=RANF(DUM) 

\ 

RNTOT=RNT0T+RN 

i 

End  Loop  ' C ' 

180  ) 

RNS . =C1 ( RNT0T=6 . 0 ) +C2 

^ \ 

End  Loop  'B' 


SUBROUTINE  'ELEM* 


Purpose: 


Method : 


To  compute  matrix  elements  in  one  of  two  ways 
specified  by  'N' 


N=0 


x(l,l) 

X(2,l) 

x(3,l) 

x(l,2) 

x(2,2) 

x(3,2) 

x(l,3) 

x(2,3) 

x(3,3) 


N=1 


x(l,l) 

x(2,l) 

x(3,l) 

x(l,2) 

x(2,2) 

z(3,2) 

x(l,3) 

x(2,3) 

x(3,3) 


cos(v(<)  cos  (9) 

-$in(4;)  cos((j))  )-  cos(i(j)  sin(9)  sin((t>) 
sinCij;)  sin((t))  + cos(i^)  sin(6)  cos((t)) 
sin(<|^)  cos(0) 

cos((i;)  cos(0)  + sin(i|;)  sin(0)  sin((})) 

-cos(4))  sin((|))  + sin(iij)  sin(0)  cos((j)) 

-sin(9) 

cos(0)  sin(())) 

cos(0)  COS((t)) 


cos(i|;)  cos  (9) 
sin(tjj)  cos(0) 

-sin (0) 

-sin(i|j)  cos(()))  + cos(\}i)  sin(9)  sin(())) 
cos(i|^)  cos((>)  + sin(ii;)  sin(0)  sin((()) 
cos(0)  sin((i)) 

sin(ti;)  sin((i))  + cos('i{))  sin(9)  cos((l>) 
-cos(4')  sin(<{))  + sin(ij;)  sin(0)  cos((|)) 
cos(0)  COS(())) 
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Purpose;  Given  the  matrices  X,  Y,  AT,  AA — Solve 
for  A: 

A = X*Y*AT-AA  Where 

X=3x3 

Y=3x3 

AT=3xl 

AA=3xl 


Method ; 


a) 

C.  = Y. ,1 
1 1 

AT,  + Y. ,2  AT-  + Y. ,3  AT- 
1 1 2 1 ' 3 

(i=l,3) 

b) 

D.  = X. ,1 
1 1 

+ X^,2  C2  + X^,3  C3 

(i=l,3) 

c) 

A.  = D.  - 
1 1 

AA^  (i=l,3) 
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APPENDIX  D 


RADAR  LAG  PROGRAM 


The  SDATA  program  requries  two  cards  of  input  for  each 
pass  which  is  to  be  analyzed.  The  first  data  card  contains 
the  following  input  which  directs  the  reading  of  the  Sight  Eval 
tape  being  used. 


FlO.l  - XMISN  - The  mission  number 

110  - IPAS  - The  pass  to  be  analyzed  from  the  above 

mission 

110  - ITME  - That  time  during  the  pass  at  which  compu- 

tations are  to  begin  (msec) 

110  - NPTS  - The  number  of  data  time  intervals  over 

which  SDATA  is  to  perform  computations 


Twenty-five  Sight  Eval  data  arrays  are  required  as  input 
to  SDATA.  Each  of  these  arrays  is  smoothed  by  SDATA.  The 
number  of  points  to  be  used  in  each  smooth  must  be  specified. 
This  is  read  in  on  the  second  data  card. 


2512  - NFIT(I) ,1=1,25 


The  number  of  points  to  be  used 
in  smoothing  each  of  the  25  Sight 
Eval  variable  arrays  input  to 
SDATA . 
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